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Quantum-confined semiconductor structures are the cornerstone of modern-day electronics. Spa- 
tial confinement in these structures leads to formation of discrete low-dimensional subbands. At 
room temperature, carriers transfer among different states due to efficient scattering with phonons, 
charged impurities, surface roughness and other electrons, so transport is scattering-limited (dif- 
fusive) and well described by the Boltzmann transport equation. In this review, we present the 
theoretical framework used for the description and simulation of diffusive electron transport in 
quasi-two-dimensional and quasi-one-dimensional semiconductor structures. Transport in silicon 
MOSFETs and nanowires is presented in detail. 
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I. INTRODUCTION 

Quantum-confined semiconductor structures, 
such as silicon MOSFETa or GaAs-based resonant- 
tunneling diodes, are the cornerstone of modern-day 
electronics Spatial confinement in these structures 
leads to formation of discrete low-dimensional subbands: 
energy levels are quantized in each direction of confine- 
ment, while the momentum remains a good (continuous) 
quantum number in the unconfined directions. If carriers 
are confined along one direction and free to move in 
the two-dimensional (2D) plane perpendicular to it, the 
structure is being referred to as a quasi-two-dimensional 
electron gas (Q2DEG). The structure is considered a 
quasi- one- dimensional electron gas (Q1DEG) if the car- 
rier motion is unbound in one dimension (ID) as a result 
of 2D-confinement. At room temperature, transport 
within each subband and transitions among subbands 
can essentially be described semiclassically, using the 
Boltzmann transport equation. Its state-of-the-art 



solution is obtained via the ensemble Monte Carlo 
technique.— At low temperatures, quantum-coherence 
effects become prominent, and transport ceases to 
be semiclassical in nature. In this regime, transport 
description is better achieved using the Wigncr-function 
formalism^ or nonequilibrium Green's functions,— >£ as 
long as the single particle picture is valid. However, 
in the remainder of this text, we will be concerned 
with the semiclassical transport picture, appropriate for 
room-temperature operation of Q2DEG and Q1DEG 
electronic structures. 

In structures such as the resonant tunneling diode^ 
or the quantum cascade laser,— the actual device opera- 
tion is based on utilization of quantum confinement and 
tunneling. In contrast, in structures such as silicon MOS- 
FETs, quantum-confinement features emerge as a result 
of miniaturization and are usually detrimental: tunneling 
through the gate oxidep^ source-to-drain tunneling and 
space-quantization effects are expected to be important 
in nanoscale MOSFETs and HEMT devices and require 
a solution of the ID Schrodinger-Poisson problem. So- 
lution of the 2D Schrodinger-Poisson problem is needed, 
for example, for describing the channel charge in narrow- 
width MOSFETs 1 ^ and nanowires i 1 ^ 3 - 

Successful scaling of MOSFETs towards shorter chan- 
nel lengths requires thinner gate oxides and higher doping 
levels to achieve high drive currents and minimized short- 
channel effects! 14 ' 15 As the oxide thickness is scaled to 
below 10 nm, quantum confinement of inversion charge 
leads to an appreciable inversion layer capacitanc o 16 ' 17 
in series with the oxide, so the total gate capacitance is 
lowered. (Further modification of the gate oxide capac- 
itance stems from the image and many-body exchange- 
correlation effects in the inversion layer as well as poly- 
silicon gate depletion^) 

The low-field electron mobility is an important quan- 
tity that determines the performance of semiconductor 
devices. Surface roughness scattering (SRS) is by far the 
most important cause of mobility degradation in conven- 
tional MOSFETs at high transverse fields. One would 
expect the SRS to be even more detrimental in silicon 
nanowires (SiNWs) than in conventional MOSFETs be- 
cause SiNWs have four Si-Si02 interfaces, as opposed to 
one such interface in conventional MOSFETs. This has 
indeed been confirmed recently on ultrathin cylindrical 12 
and rectangular nanowiresi^ 3 - In addition, confinement 
leads to a modification in the acoustic phonon spectrum 
in SiNWs ; 13 ' 20 which leads to increased electron-phonon 
scattering and a lowered electron mobilit y 21 ' 22 ' 23 

In this review, we present the theoretical framework 
typically used for the description and simulation of dif- 
fusive electron transport in quasi-2D and quasi-lD semi- 
conductor structures. In Sec. UH we overview the for- 
mation of the Q2DEG in Si inversion layers and GaAs- 
modulation doped heterostructures. We discuss the so- 
lution to the ID Schrodinger-Poisson problem in the di- 
rection of confinement and the density of states calcula- 
tion. In Sec. IIIIl we overview the scattering mechanisms 



in semiconductors, their origin and model Hamiltonians, 
and give the scattering rates for Q2DEGs in the Born ap- 
proximation. In Sec. IIV1 we present the ensemble Monte 
Carlo simulation of the electron mobility in the inver- 
sion layer Q2DEG formed near the Si/Si02 interface of a 
nanoscale MOSFET. Section[V]introduces the Q1DEG in 
nanowires, and overviews phonon confinement and band- 
structure modification in these structures. In Sec. IVI1 
scattering rates for the Q1DEG are given, while the elec- 
tron mobility in thin silicon nanowires, as obtained from 
detailed ensemble Monte Carlo simulations, is presented 
in Sec. EE3] 



II. QUASI-2D ELECTRON SYSTEMS 
A. Silicon Inversion Layers 

The best known examples of the Q2DEG are silicon 
MOSFETs and GaAs/AlGaAs heterostructures. An in- 
tegral part of any MOSFET device is the metal-oxide- 
semiconductor (MOS) capacitor. Regarding the MOS 
capacitors, the induced interface charge is closely linked 
to the shape of the electron energy bands of the semi- 
conductor near the interface. At zero applied voltage, 
the bending of the energy bands is ideally determined by 
the difference in the work functions of the metal and the 
semiconductor. This band bending changes with the ap- 
plied bias and the bands become flat when we apply the 
so-called flat-band voltage given by 

VfB = $M - (1) 

where $Af and Q sc are the work functions of the metal 
and the semiconductor, respectively. The various ener- 
gies involved are indicated in Fig. [TJ where we show 
typical band diagrams of an MOS capacitor at zero bias. 
Xs is the electron affinity for the semiconductor, E c is 
the energy of the conduction band edge, and Ep is the 
Fermi level at zero applied voltage. 

In stationary conditions, no net current flows in the 
direction perpendicular to the interface, owing to the 
very high resistance of the insulator layer (however, this 
does not apply to very thin oxides of a few nanometers, 
where tunneling becomes important). Hence, the Fermi 
level will remain constant inside the semiconductor, irre- 
spective of the biasing conditions. However, between the 
semiconductor and the metal contact, the Fermi level is 
shifted by Epu — Ep s = €Vq (see Figs. [2]and[3J). Hence, 
we have a quasi-equilibrium situation in which the semi- 
conductor can be treated as if in thermal equilibrium. 

An MOS structure with a p-type semiconductor will 
enter the accumulation regime of operation when the 
voltage applied between the metal and the semiconduc- 
tor is more negative than the flat-band voltage (Vfb < 
in Fig. [5]). In the opposite case, when Vg > Vfb, the 
semiconductor-oxide interface first becomes depleted of 



2 



Ideal MOS capacitor with Real MOS capacitor with 

a p-type semiconductor a .p-type semiconductor 




FIG. 1. (Left panel) Energy band diagram of an ideal MOS 
capacitor. (Right panel) Energy band diagram of a real MOS 
capacitor. 
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FIG. 2. MOS capacitor under accumulation conditions. 



holes and we enter the so-called depletion regime (Fig. [3J 
top panel). By increasing the applied voltage, the band 
bending becomes so large that the energy difference be- 
tween the Fermi level and the bottom of the conduction 
band at the insulator-semiconductor interface becomes 
smaller than that between the Fermi level and the top of 
the valence band. This is the inversion regime of opera- 
tion (Fig. [31 bottom panel). 

Carrier statistics tells us that the electron concentra- 
tion will then exceed the hole concentration near the in- 
terface and we enter the inversion regime. At a larger 
still applied voltage, we finally arrive at a situation in 
which the electron density at the interface exceeds the 
doping density in the semiconductor. This is the strong 
inversion case, in which we have a significant conducting 
sheet of inversion charge at the interface (Fig. [31 bottom 
panel). In the description that follows, symbol <p is used 
to denote the potential in the semiconductor measured 
relative to the potential at a position x deep inside the 
semiconductor. Note that tp becomes positive when the 




FIG. 3. MOS capacitor under depletion conditions (top 
panel) and inversion conditions (bottom panel). 

bands bend down, as in the example of a p-type semi- 
conductor shown in Figure [3j From equilibrium electron 
statistics, we find that the intrinsic Fermi level Ei in the 
bulk corresponds to an energy separation etpF from the 
actual Fermi level Ep of the doped semiconductor, 

i p F = V T ln(^\ >0, (2) 

where Vt is the thermal voltage, Na is the shallow ac- 
ceptor density in the p-type semiconductor and rii is 
the intrinsic carrier density of silicon. According to the 
usual definition, strong inversion is reached when the to- 
tal band bending equals 2etp p , corresponding to the sur- 
face potential (p s = 2ipp . Values of the surface potential 
such that < (p s < 2ifF correspond to the depletion and 
the weak inversion regimes, respectively, tp s — is the 
flat-band condition, and ip a < corresponds to the accu- 
mulation mode. Note that, deep inside the semiconduc- 
tor, we have y(oo) = 0. Under the flat-band condition 
(Vg — Vfb): the surface charge is equal to zero. In ac- 
cumulation (Vg < Vfb)-, the surface charge is positive, 
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and in depletion and inversion Vg > Vfb), the surface 
charge is negative. 

In order to relate the semiconductor surface poten- 
tial <p s to the applied voltage Vg, we have to investigate 
how this voltage is divided between the insulator and the 
semiconductor. Using the condition of continuity of the 
electric flux density at the semiconductor-insulator inter- 
face, we find e s F s — e ox F ox , where e ox and e s are the ab- 
solute permittivities (also known as dielectric constants) 
of the oxide layer and the semiconductor, respectively, 
while F ox and F s are the respective electric fields in the 
two materials. Hence, for an insulator of thickness d ox , 
the voltage drop across the insulator becomes F ox d ox . 
Accounting for the flat-band voltage, the applied voltage 
can be written as Vg — Vfb + fs + £sF s /C ox , where 
C ox = £ox/d ox is the insulator capacitance per unit area. 

The threshold voltage, Vth, is the gate voltage cor- 
responding to the onset of strong inversion. It is one 
of the most important parameters characterizing metal- 
insulator-semiconductor devices. As discussed above, 
strong inversion occurs when the surface potential, tp s , 
becomes equal to 2ipp- For this surface potential, the 
free charge induced at the insulator-semiconductor inter- 
face is still small compared to the charge in the depletion 
layer, and the threshold voltage is calculated using: 



Vg 



tp s H — , where C 



Vth — Vg for which ip s = 2(pp. 



(3) 



Note that the threshold voltage may also be affected by 
the so-called fast surface states at the semiconductor- 
oxide interface and by fixed charges in the insulator layer. 
However, this is not a significant concern with modern- 
day fabrication technology. 

The threshold voltage separates the subthreshold 
regime, where the mobile carrier charge increases expo- 
nentially with increasing applied voltage Vg, from the 
above-threshold regime, where the mobile carrier charge 
is linearly dependent on the applied voltage Vg- How- 
ever, there is no clear point of transition between the two 
regimes, so different definitions and experimental tech- 
niques have been used to determine Vth- Well above 
threshold, the charge density of the mobile carriers in 
the inversion layer can be calculated using the parallel- 
plate charge control model. This model gives an ade- 
quate description for the strong inversion regime of the 
MOS capacitor, but fails for applied voltages near and 
below threshold (i.e., in the weak inversion and deple- 
tion regimes). Several expressions have been proposed 
for a unified charge control model (UCCM) that covers 
all regimes of operation. 

Successful scaling of MOSFETs towards shorter chan- 
nel lengths requires thinner gate oxides and higher dop- 
ing levels to achieve high drive currents and minimized 
short-channel effects! 14 ' 15 For these nanometer devices it 
was demonstrated that, as the oxide thickness is scaled 
to 10 nm and below, the total gate capacitance is lower 
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FIG. 4. SCHRED 30 simulation data for the shift in the 
threshold voltage compared to the experimental values pro- 
vided by van Dort and co-workers . 27 ' 28 



than the oxide capacitance due to the comparable val- 
ues of the oxide and the inversion layer capacitances. As 
a consequence, the device transconductance is degraded 
relative to the expectations of the scaling theory^ 

The two physical origins of the inversion layer capac- 
itance, the finite density of states and the finite inver- 
sion layer thickness, were demonstrated experimentally 
by Takagi and Toriumi. 16 A computationally efficient 
three-subband model, that predicts both the quantum- 
mechanical effects in the electron inversion layers and 
the electron distribution within the inversion layer, was 
proposed and implemented into the PISCES simulator. 17 
The influence of the image and many-body exchange- 
correlation effects on the inversion layer and the total 
gate capacitance was studied by Vasileska et al*& It was 
also pointed out that the depletion of the poly-silicon 
gates considerably affects the magnitude of the total gate 
capacitance ^ 

The inversion layer capacitance was also identified as 
being the main cause of the second-order thickness de- 
pendence of the MOSFET IV-characteristicsi^ 5 . The finite 
inversion layer thickness was estimated experimentally 
by Hartstein and Albert The high levels of substrate 
doping were found responsible for the increased threshold 
voltage and decreased channel mobility. A simple analyt- 
ical model that accounts for this effect was proposed by 
van Dort and co-worker o 27 ' 28 and confirmed by Vasileska 
and Ferry 2 ^ by investigating the doping dependence of 
the threshold voltage in MOS capacitors (Fig. 0J. 



B. ID Schrodinger-Poisson Problem for Silicon 
Inversion Layers 

The periodic crystal potential in bulk semiconducting 
materials is such that, for a given energy in the conduc- 
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tion band, the allowed electron wavevectors trace out a 
surface in k-space. In the effective-mass approximation 
for silicon, these constant energy surfaces can be visu- 
alized as six equivalent ellipsoids of revolution (Fig. O, 
whose major and minor axes are inversely proportional 
to the effective masses. A collection of such ellipsoids for 
different energies is referred to as a valley. 

In this framework, the bulk Hamiltonian for an elec- 
tron, residing in one of these valleys is of the form 

( ^ ^ ^ ^ ^ ^ \ t r / \ 

Ho{R) -{^ x d^ + 2^ y d^ + 2^ z d^) +Vcs{z) 
= H 0{l (r) + H ox (z), (4) 

where R = (r, z), V e s(z) — Vh{z) + V cxc {z) is the ef- 
fective potential energy profile of the confining poten- 
tial along the z-direction, Vh(z) is the Hartree potential 
which is nothing more but a solution of the ID Pois- 
son equation introduced later in the text, V cxc (z) is the 
exchange-correlation potential (also discussed later in the 
text), i? || is the parallel part of Ho (associated with the 
motion in the xy-plane, perpendicular to the confinement 
direction), and the transverse part is defined as 

^(z) = -^^ + V Bff (z). (5) 



The basis-states of the unperturbed Hamiltonian are as- 
sumed to be of the form 

V„(R) = 4f e ik ' r Vn(2), (6) 

where k is a wavevector in the rry-plane and A is the 
area of the sample interface. The subband wavefunctions 
satisfy the one-dimensional Schrodinger equation 

H 0± (z)ip n (z) = £ n ip n (z) (7) 

subject to the boundary conditions that tp n {z) are zero 
for z — and approach zero as z — > 00. In Eq. (|7|). 
£ n is the subband energy and ?p n (z) is the corresponding 
wavefunction. In the parabolic band approximation, the 
total energy of an electron is given by 

h 2 k 2 

£ n (k) = — r +£ n =£ k + £ n , (8) 



where £k is the kinetic energy and m* y is the density-of- 
states mass along the xy-plane. 



An accurate description of the charge in the inversion 
layer of deep-submicrometer devices and, therefore, the 
magnitude of the total gate capacitance, C'tot, requires a 
self-consistent solution of the ID Poisson equation 



d_ 

dz 



e(z) 



dip 



[N+(z)-NX(z)+p(z)-n(z)] 



and the ID Schrodinger equation 



(9) 



if 



2mjr dz 2 



+ V eS (z) 



ip V j(z) = £vj1pv 3 {z). (10) 



In Eqs. ([9]) and Ijl0|> . ip(z) is the electrostatic poten- 
tial [the Hartree potential Vh(z) = — e(f(z)], e(z) is the 
spatially dependent dielectric constant, Np and N^(z) 
are the ionized donor and acceptor concentrations, n{z) 
and p(z) are the electron and hole densities, V e f{ is the 
effective potential energy term that equals the sum of 
the Hartree and exchange-correlation corrections to the 
ground state energy of the system, is the effective 
mass normal to the semiconductor-oxide interface of the 
f-th valley, while £ V j and ip V j(z) are the energy level and 
the corresponding wavefunction of the electrons residing 
in the j-th subband from the w-th valley. The electron- 
density is calculated using 



n(z) 



N vj \ip v j{z)\' 



(11) 



where N v j is the sheet electron concentration in the j-th 
subband from the f-th valley, given by 
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N vj = Vv ^k B T\n{l + exp [(E F - £ VJ )/k B T}} (12) 



where v v is the valley degeneracy factor and Ep is the 
Fermi energy. When evaluating the exchange-correlation 
corrections to the chemical potential, we have relied on 
the validity of the density functional theory (DFT) of 
Hohenberg and Kohn^ and Kohn and Sham^ Accord- 
ing to DFT, the effects of exchange and correlation can 
be included through a one-particle exchange-correlation 
term Vcxc [n(z)], defined as a functional derivative of 
the exchange-correlation part of the ground-state en- 
ergy of the system with respect to the electron den- 
sity n(z). In the local density approximation (LDA), 
one replaces the functional V^ xc [n(z)] with a function 
V cxc [n(z)] = fi cxc [no =n(z)], where (i cxc is the exchange- 
correlation contribution to the chemical potential of a 
homogeneous electron gas of density uq, which is taken 
to be equal to the local electron density n(z) of the in- 
homogeneous system. In our model, we use the LDA 
and approximate the exchange-correlation potential en- 
ergy term 14 xc [n(z)] by an interpolation formula devel- 
oped by Hedin and Lundqvisti^ Exchange and correla- 
tion effects tend to lower the total energy of the system, 
and, as discussed later, lead to a non-uniform shift of the 
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FIG. 5. (Left panel) Constant-energy surfaces for the conduction-band of silicon, showing six conduction-band valleys in the 
(100) direction of momentum space. The band minima, corresponding to the centers of the ellipsoids, are 85% of the way to 
the Brillouin-zone boundaries. The long axis of an ellipsoid corresponds to the longitudinal effective mass of the electrons in 
silicon, mi = 0.916mo, while the short axes correspond to the transverse effective mass, rat = 0.190mo (mo is the free-electron 
mass). For a (100) surface, the A2-band has the longitudinal mass (mi) perpendicular to the semiconductor interface and the 
A4-band has the transverse mass (m t ) perpendicular to the interface. Since larger mass leads to a smaller kinetic term in the 
Schrodinger equation, the unprimed lader of subbands (as it is usually called), corresponding to the A2-band, has the lowest 
ground state energy. The degeneracy of the unprimed ladder of subbands for a (100) surface is 2. For the same reason, the 
ground state of the primed ladder of subbands corresponding to the A4-band is higher than the lowest subband of the unprimed 
ladder of subbands. The degeneracy of the primed ladder of subbands for a (100) surface is 4. (Right panel) Potential diagram 
for inversion of p-type silicon. The notation vj refers to the j-ih subband of either the A2-band (v — 1) or A-4-band (v = 2). 



energy levels and repopulation of the various subbands. 
The enhancement of the exchange-correlation contribu- 
tion to the energy predominantly affects the ground sub- 
band of the occupied valley; the unoccupied subbands of 
the same valley are essentially unaffected. As a result, a 
noticeable increase in the energy of intersubband transi- 
tions can be observed at high electron densities (more on 
electron scattering can be found in Sec. 



C. GaAs/AlGaAs Heterostructures. Effective 
Mass Schrodinger Equation for Heterostructures 

While our primary focus in this paper will remain on 
Si-based electron system, in this section we will briefly 
discuss another important Q2DEG: the modulation- 
doped GaAs-AlGaAs heterostructures The bandgap in 
AlGaAs is wider than in GaAs. By variation of doping 
it is possible to move the Fermi level inside the forbid- 
den gap. When the materials are put together, a unified 
chemical potential is established, and an inversion layer 
is formed at the interface (see Fig. [5]). 



The Q2DEG created by modulation doping can be 
squeezed into narrow channels by selective depletion in 
spatially separated regions. The simplest lateral confine- 
ment technique is to create split metallic gates in a way 
shown in Fig. [7] We refer the reader to extensive litera- 
ture on transport properties of GaAs/ AlGaAs and other 
modulation doped heterostructures (for instance, see Ref. 
1- 

In semiconductors, some of the most interesting appli- 
cations of the Schrodinger equation in the effective mass 
approximation involve spatially varying material compo- 
sitions and heterojunctions. The effective mass approx- 
imation can still be used with some caution. Since the 
effective mass is a property of the bulk, it is not well de- 
fined in the neighborhood of a sharp material transition. 
In the hypothesis of slow material composition variations 
in space, one can adopt the Schrodinger equation with a 
spatially varying effective mass, taken to be the mass of 
the bulk with the local material properties. However, it 
can be shown that the Hamiltonian operator is no longer 
Hermitian for varying mass. A widely used Hermitian 
form brings the effective mass inside the differential op- 
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erator as — ^-V • (-^rV^>). This approach is extended to 
abrupt heterojunctions, as long as the materials on the 
two sides have similar properties and bandstructure, as 
in the case of the GaAs/AlGaAs system with less than 
45 % Al. One has to keep in mind that very close to the 
heterojunctions the effective mass Schrodinger equation 
provides a reasonable mathematical connection between 
the two regions, but the physical quantities arc not neces- 
sarily well defined. For instance, in the case of a narrow 
potential barrier obtained by using a thin layer of Al- 
GaAs surrounded by GaAs, it is not clear at all what 
effective mass should be used for the AlGaAs, since such 
a region can certainly not be approximated by the bulk. 
It is even more difficult to treat the case when there is a 
transition between direct and indirect bandgap materials 
(for example, GaAs and AlGaAs with over 45% Al). 

Assuming a uniform mesh size Ax, the Hamiltonian 
of the Schrodinger equation can be discretized in ID by 
introducing midpoints in the mesh on both sides of a 
generic grid point i. First, we evaluate the outer deriva- 
tive at point i with centered finite differences, using quan- 
tities defined at points (i — 1/2) and (i + 1/2) 
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FIG. 6. Band diagrams near the interface between n-AlGaAs 
and intrinsic GaAs, (a) before and (b) after the charge trans- 
fer. 
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FIG. 7. On the formation of a narrow channel by a split gate. 
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and then the derivatives defined on the midpoints are 
also evaluated with centered differences using quantities 
on the grid points 



ti 2 



2(Axf 



ip(i + 1) - ip(i) ip(i) - ip(i - I) 
m*(i + l/2) m*(i - 1/2) 



• (14) 



The effective mass is the only quantity which must be 
known at the midpoints. If an abrupt heterojunction is 
located at point i, the abrupt change in the effective mass 
is treated without ambiguity. 



D. Density of States (DOS) for Low-Dimensional 
Systems 

An important quantity characterizing a quantum- 
mechanical system is the density of states (DOS) func- 
tion. The density of states g{£) is defined as the number 
of states per energy interval (£, £ + d£). It is clear that 



g(£) = J2S(£-£ a ), 



(15) 



where a is the set of quantum numbers characterizing the 
states. In the present case, it includes the subband quan- 
tum number n, spin quantum number a, valley quantum 
number v and the in-plane quasi-momentum k. If the 
spectrum is degenerate with respect to spin and valleys, 



7 



one can define the spin degeneracy v s and the valley de- 
generacy v v to get 



v s v v \ - 



(16) 



Here we calculate the number of states per unit volume, 
D being the dimensionality of the space. For a 2D case, 
we obtain 



52e {£-£*). 



(17) 



Within a given subband, the 2D density of states func- 
tion is energy-independent. Since there can exist several 
subbands in the confining potential, the total density of 
states can be represented as a set of steps, as shown in 
Fig. [5] At a low temperature (fc^T <C Ep), all the states 
are filled up to the Fermi level. Because of the energy- 
independent density of states, the sheet electron density 
is linear in the Fermi energy, namely 



N„ = N 



v s v v m* Ep 
(2Tr)h 2 



(18) 



The Fermi momentum in each subband can be deter- 
mined as 



y/2m(E F - £ n ) 



(19) 



In Eq. (|18p , N is the number of transverse modes having 
the edge £ n below the Fermi energy. 

The situation is more complicated if the gas is further 
confined in a narrow channel, say, along the y-axis. The 
in-plane wave function can be decoupled as a product 



the corresponding energy being 



2 1,2 



£n ~b £s 



2m 



(20) 



(21) 



In the last equation, £ ns = £ n + £ s characterizes the 
energy level in the potential confined in both (z and y) 
directions. For square-box confinement, the terms are 



£ a 



(sirh) 2 
2mW 21 



(22) 



where W is the channel width, while for parabolic con- 
finement U (y) 
structures), we have 



(l/2)mcu 2 y 2 (typical for split-gate 



£ s = (s — 1/2) Tiujq 



(23) 



For these systems, confined in 2D, the total density of 
states is 



g{£) 



2 3 / 2 nh 



E 



e(£ - £ns) 



(24) 



The energy dependence of the density of states is shown 
in Fig. M 
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FIG. 8. Density of states for a quasi-2D system. The three 
lowest subbands are included. 
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FIG. 9. Density of states for the bulk (3D, blue), quantum 
well (2D, red), quantum wire (ID, green) and quantum dot 
(0D, black). 



III. SCATTERING IN QUASI-2D ELECTRON 
SYSTEMS 



Charge transport in the diffusive regime is governed 
by carrier scattering from lattice vibrations, charged im- 
purities, defects, interface roughness, as well as other 
electrons. Calculation of the scattering rates for con- 
fined carriers proceeds in a similar manner as in the 3D 
casej^iM but proper wavefunctions for 2D carriers must 
be used. Before we go into the details of the calculation 
of the matrix elements of some of the most important 
scattering mechanisms listed in Fig. I10I we will derive a 
few expressions. 
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FIG. 10. Scattering mechanisms in a typical semiconductor. 



Suppose we want to calculate the scattering rate out 
of some state k in subband n. For that purpose we will 
use Fermi's golden rule, which gives us the transition 
rate from state k in subband n into state k' belonging to 
subband m by means of emission of an energetic particle 
(e.g. a phonon or a photon) with energy Tiuj : 

S„ m (k,k') = — \M(k,-k')\ 2 nm 5(£' -£ + hu). (25) 

Assuming a plane- wave basis for the wavefunctions in the 
unconfined direction (xy-plane), the total wavefunctions 
of the initial and the final states are of the following gen- 
eral form (for a Q2DEG): 



V„(k, z) = 4=e 4k '>„(z) , V™(k', z') = 4=e lk ' 
VA VA 



(26) 

where A is the area of the sample, r is the position vector 
in the xy-plane and R = (r, z) is a 3D position vector. 
The matrix element for scattering between states k and 
k in subbands n and m, respectively, is then given by 

M(k,k')„ m = I Je« k - k '^d 2 r J y* m {z)H qv (Tl)y n {z) 

where H qv is the interaction potential and the form of the 
integral with respect to z depends upon the type of the 
scattering dynamics consideredi 2 ' 34 i 35 In low-dimensional 
systems, since the momentum is quantized in one ore two 
directions to form subbands, it is important to note that 
we now have additional intrasubband and intersubband 
transitions, which significantly complicates the genera- 
tion of the scattering tables and choosing the final state 
after scattering. Below, we give the matrix elements for 
some of the most important scattering mechanisms that 



are present in silicon inversion layers and GaAs/AlGaAs 
heterostructures. 



A. Electron-phonon scattering 

Phonon scattering can cause three different types of 
electronic transitions in the Si inversion layer: transi- 
tions between states within a single valley via acous- 
tic phonons (called intravalley acoustic-phonon scat- 
tering) and nonpolar optical phonons (called intraval- 
ley optical phonon scattering), and transitions between 
different valleys mediated by high-momentum acoustic 
or nonpolar optical phonons (called together interval- 
ley scattering) ^L^^^^. Intravalley acoustic- 
phonon scattering involves phonons with low energies 
and is an almost elastic process. The intravalley optical- 
phonon scattering is induced by optical phonons of low 
momentum and high energy. Intervalley scattering can 
be induced by the emission and absorption of high- 
momentum, high-energy phonons, which can be of cither 
acoustic- or optical-mode variety. Intervalley scattering 
can therefore be important only for temperatures high 
enough that an appreciable number of suitable phonons 
is excited or for hot electrons that can emit high energy 
phonons^ l n order to evaluate the scattering potential 
that describes the electron-phonon interaction, we need a 
Hamiltonian that describes the coupled electron-phonon 
system. The total Hamiltonian of the system is given 
byi^ 



H — H e + H a + H e 



(28) 



where H e is the electronic part, H a is the atomic part 
that describes the normal modes of vibration of the solid, 
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H, 



= J2 V *» ( r *- R ;)- 



(29) 



In general, each ion is at a position Rj = R^+Uj, which 

is a sum of the equilibrium position R^ and the displace- 
ment Uj. Under the assumption of small displacements, 
one can expand V ea in a Taylor series 

V ea (r-Rf) = Veair-TLf^-UyVrVeair-nf ) + o(Q 2 ). 

( 30 ) 

The zero-order term is the potential function for the elec- 
trons when the atoms are in their equilibrium positions, 
which forms a periodic potential in the crystal. The so- 
lution of the Hamiltonian for electron motion in this pe- 
riodic potential gives the Bloch states of the solid. Since 
the first order term is much smaller than the zero-order 
term, the electron-phonon interaction can be treated per- 
turbatively. Therefore, the lowest order term for the 
electron-phonon interaction is of the form 



where MN = pV, and p is the density of the solid. The 
summation over G represents summation over all recip- 
rocal lattice vectors of the solid. If one defines a function 



M, 



qA 



(— 



5> 4G - r (q + G).e qA K Q (q + G), 



G 



(35) 

then the Hamiltonian for the electron-phonon interaction 
becomes 



H e - ph (r) 



M qA e 



iq-r 



(36) 



The exact form of the matrix elements for acoustic 
and nonpolar-optical phonon scattering (zero- and first- 
order terms) are given below. Since we will need to 
make a clear distinction between 3D and 2D vectors, 
in what follows we are going to use the following nota- 
tion: capital bold letters will refer to three-dimensional 
vectors, whereas small bold letters will be used for two- 
dimensional wavevectors that lie in the xy-plane. 



J2 u r V r T4 a (r-R} 0) ) 



(31) 



It is obvious that this interaction Hamiltonian does not 
act on the spin variables in this approximation. The 
Fourier transform of V ea can be written as 



V ea (r) 



(32) 



where N is the number of primitive cells, and wavevector 
q spans the whole q-space. Ionic displacement may be 
decomposed into normal-mode representation and it is 
customary to write 



E 



\2MNu) k \ 



1/2 



(33) 



, akAe k c 



k-R ' 



it 



l kA e k 



e^c 



ik-RV 



where ekA is the unit polarization vector that obeys the 
standard orthonormality and completeness relations, WkA 
is the phonon frequency of phonon branch A for wavevec- 
tor k (running over the whole Brillouin zone), fikA ( a k\) 
are the phonon annihilation (creation) operators. In 
acoustic waves, Uj refers to the relative displacement of 
the unit cell as a whole with respect to adjacent unit cells; 
in optical waves it refers to the relative displacement of 
the basis atoms within the unit cell. Thus 

H e - P h(r) = ^e l ^ +G )- r U eQ (q + G) (34) 

h 



x (q + G) • e qA 



a q A + a qA 



1. Deformation potential scattering 

In general, the application of mechanical stress al- 
ters the band structure by shifting energies, and, where 
it destroys symmetry, by removing degeneracies. It is 
usually assumed that the mechanical stress does not 
change the band curvature, and therefore does not change 
the effective masses, but introduces ashift in the energy 
states that are close to the band extremum i 47 i 48 i 49 i 50 
For isotropic elastic continuum, the matrix element for 
deformation potential scattering (acoustic phonons) can 
be obtained by taking the long-wavelength limit of Eq. 
([55)1 . 51 For small values of Q, the summation over recip- 
rocal lattice vectors can be neglected, except for the term 
G = 0. The screened electron-ion interaction becomes a 
constant which is usually denoted as S (it gives the shift 
of the band edge per unit elastic strain). Under these 
assumptions, Mq\ simplifies to 



Mqx = iQ ■ e QA 



/ HE 2 



\ 2pVuoQX 



(37) 



where wqa = V 8 \Q, v s \ is the sound velocity, is the 
phonon frequency. Long wavelength acoustic phonons 
(LA mode) have Q||eQ A which makes the matrix ele- 
ment non-zero. TA phonons have Q_LeQ A which makes 
the matrix element vanish. Therefore, the deformation 
potential mainly couples electrons to LA phonons. 

For anisotropic elastic continuum such as silicon, the 
deformation potential constant S becomes a tensor. The 
anisotropy of the intravalley deformation potential in 
the ellipsoidal valleys in silicon has been extensively 
studied by Herring and Vogt. 50 Expanding the electron- 
phonon matrix element over spherical harmonics and re- 
taining only the leading terms, they have expressed the 



10 



\y 




FIG. 11. Angular dependence of the deformation potential 
for longitudinal modes. 

anisotropy of the interaction in terms of the angle 9q be- 
tween wavevector Q of the emitted (absorbed) phonon 
and the longitudinal axis of the valley. They have shown 
that the matrix element is proportional to Q via the 
deformation potential A\(8q) (A =LA or TA) given 

u.,52.53.54 



A LA (6 Q )^Z d + Z u cos 2 (6 Q ), (38) 

and 



A TA (9 Q ) w E u cos(0 Q ) sin(0 Q ). (39) 

Equation (|3"9")) accounts for the contribution of both TA 
branches. Therefore, the acoustic mode scattering is 
characterized by two constants: S u (uniaxial shear po- 
tential) and 3^ (dilatation potential) that is believed to 
have values of approximately 9.0 eV and -11.7 eV, respec- 
tively. In bulk silicon, this anisotropy is usually ignored 
by using an effective deformation potential constant 3f^ 
for the interaction with longitudinal modes, and ignoring 
the role of the lower-energy TA modes. This approxima- 
tion can be justified due to the following reasons: The 
acoustic modes are most effective at low energy. In this 
regime and in the usual elastic and equipartition approx- 
imation (described later), due to the linear dependence 
on Q, scattering of electrons at some energy E samples al- 
most uniformly the constant energy ellipsoid, so that one 
can take the average values of A\ over the ellipsoid. Since 
there is nothing to fix the energy scale in the problem, 
this averaging procedure is independent of the electron 
energy. Moving to the two-dimensional situation, one 
cannot follow a parallel path to arrive at an isotropic, en- 
ergy independent effective deformation potential, which 
complicates the treatment of this scattering process. 

Since the wavefunctions of the initial and final states 
are usually expressed as a product of a one-electron Bloch 
wavefunction and a harmonic oscillator wavefunction, af- 
ter the averaging over the phonon states is performed, 



the terms inside the brackets of Eq. (|36|) that repre- 
sent phonon absorption (term dq\) and phonon emission 

(term a* x ) processes reduce to \/Nq\ and ^Nq\ + 1, 
respectively. In thermal equilibrium with a lattice at 
temperature T, the phonon occupation number Nq\ is 
given by the Bose-Einstein statistics 

N ^ = c t^Lt_ V ( 4 °) 

where Jcb is the Boltzmann constant. At high enough 
temperatures, the acoustic phonon energies are much 
smaller than the thermal energy of electrons. Therefore, 
one can expand the exponent in the denominator of Eq. 
(|40p into a series and, in the equipartition approximation, 
appropriate at high temperatures, we have 

7 rri 

Nqx + 1 ~ N Q x « » 1. (41) 

nujQ\ 

Incorporating these terms as well as the exponential term 
e iQzZ m t th e definition of Mq\, after a straightforward 
calculation one finds that the matrix element squared for 
scattering between subbands n and m due to acoustic 
phonons (for both absorption and emission processes to- 
gether), and after the averaging over Q z is performed, 
reduces to 

|(n|[/r(q)M| 2 = -^[A e A ff „J 2 F nm , (42) 
where 

POO 

F nm = / dz^ n {z)^ m {z). (43) 
Jo 

The effective deformation potential constant is calcu- 
lated from 

[A^J 2 = 7^/ dq z Al(d Q )\F nm (Q z )\ 2 , (44) 

-Trim JQ 

where 

poo 

F nm {Q z )= dzil) n {z)i) m {z)e iQ * z . (45) 
Jo 

The form- factor F nm (q z ) introduces an energy scale in 
the problem by fixing the fuzzy component, the wavevec- 
tor . This result is an expected one and follows immedi- 
ately from the uncertainty principle AzAp z > h/2. Since 
the electrons are frozen into their wavefunctions, and can- 
not oscillate in the quantized direction, the uncertainty in 
the particle's location along the z-axis has been reduced. 
Therefore, there must be a corresponding increase in the 
uncertainty in the particle's ^-directed momentum^ 
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2. Nonpolar optical phonon scattering 

The scattering of electrons by zone-center optical and 
intervalley phonons in semiconductor crystals has been 
treated rather extensively by Ferry . 44 ' 55 ! 56 The nonpolar 
optical interaction is important for intrasubband scat- 
tering as well as for scattering of electrons (and holes) 
between different minima of the conduction (or valence) 
band. The latter interaction is important for scatter- 
ing of carriers in semiconductors with many- valley band 
structure, such as Si and Ge, and in the Gunn effect, 
where scattering occurs between different sets of equiva- 
lent minima. Harrison 49 pointed out that the nonpolar 
optical matrix element may be either of zero or higher or- 
der in the phonon wavevector. In subsequent treatments 
of electron transport in which the nonpolar interaction is 
important, only the zero-order term was considered, gen- 
erally owing to the impression that the higher order terms 
are much smaller. Although this is usually the case, there 
arise many cases in which the zero-order term is forbid- 
den by the symmetry of the state involved. In these cases, 
the first order term becomes the leading term, and can 
become significant in many instances. For example, the 
first-order intervalley scattering plays an important role 
in hot-electron transport in the n-type inversion layer 
in Si. Ignoring this scattering process means that there 
will be no saturation of the drift velocity at high electric 
fields, because the zero-order intervalley scattering rate is 
weakly dependent on the electron energy of high-energy 
electrons, while the first-order intervalley scattering rate 
increases as the electron energy increases. The matrix 
element for nonpolar optical phonon scattering is gen- 
erally found from a deformable ion model explained in 
the introduction part of this section. If one thinks of 
an optical phonon as occurring at finite G, then the Q 
dependence is unimportant, so that the entire matrix el- 
ement becomes constant and we have 



M QA = Max = 



hDl 



1/2 



(46) 



where D\ is the deformation field (usually given in 
eV/cm) and ojoa is the frequency of the relevant phonon 
mode which is usually taken to be independent of the 
phonon wave- vector for optical and intervalley processes. 
Fourier transforming back to real space, a constant in 
Q-space produces a delta-function in real space. There- 
fore, this zero-order term represents a short-ranged inter- 
action. A local dilatation or compression of the lattice 
produces a local fluctuation in the energy of the electron 
or hole. Incorporating the exponential term cxp(iQ z z) 
into the definition of Mo a, after averaging over Q z we 
find 



(n\U 



op(0) 



m) 



2pVuj x 



(47) 





a = 1 


a = 2 


a = 1 


v\ = 1; v% = 


v\ — 0; v-2 — 4 


a = 2 


vi = 0; v 2 = 2 


v\ = 1; u 2 = 2 



TABLE I. Degeneracy factors for transition between un- 
primed and primed subbands, for both g- and f-phonons. 



for the squared matrix element for scattering between 
subbands n and m that belong to the a and (3 valleys, 
respectively. When the zero-order matrix element for the 
optical or intervalley interaction vanishes, D\ is identi- 
cally zero. In this case, one has to consider the first-order 
term of the interaction whose matrix element is 



Mq\ = iQ ■ e Q A 



( &gj \ 
\2pVuox) 



1/2 



(48) 



In this context, a first order process means a process sim- 
ilar to acoustic phonon scattering. Following the previ- 
ously explained procedure, we find that the matrix el- 
ement squared for scattering between subbands n (a- 
valley) and m (/?- valley) is given by 



[n\uf l) \m) 



2pVuj \ 



(,Q F nm -\- c nm ), (49) 



where 



dz 



(50) 



The constant term c nm is a small correction term^ 5 . 

In the scattering among equivalent valleys, there are 
two types of phonons that might be involved in the pro- 
cess (see Fig. [T^]). The first type, the so-called g-phonon 
couples the two valleys along opposite ends of the same 
axis, i.e. [100] to [100]. This is an umklapp process and 
has a net phonon wavevector 0.37r/a. The f-phonons cou- 
ple a (100) valley with (010), (001), etc. The reciprocal 
lattice vector involved in the g-process is Gioo and that 
for an f-process is Gm. Degeneracy factors (g r ) for tran- 
sition between unprimed (a = 1) and primed (a = 2) set 
of subbands, for both g- (r = 1) and f-phonons [r = 2) 
are summarized in Table [T] 

Within a three-subband approximation (subband £q 
in a [100] valley and £\ in [100], and a generic subband 
£q in one of the other four valleys - (010) (001)). Scat- 
tering between the £q and E\ subbands in the two val- 
leys along the same axis involves only g-type phonons. 
Scattering between these two minima is usually treated 
by using a high-energy phonon of 750 K activation tem- 
perature (treated as zero-order interaction) and 134-K 
phonon treated via first order interaction. Scattering 
between £q or £\ and the four £q subbands involves f- 
phonons with activation temperatures of 630 K and 230 
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FIG. 12. Diagrammatic representation of intervalley transi- 
tions due to g- and f-phonons. 

K, treated via zero-order and first-order interaction, re- 
spectively. Scattering among subbands £q involves both 
g- and f-phonons with activation temperatures of 630 K 
(zero-order interaction) and 190 K (first-order interac- 
tion). All of the high energy phonons are assumed to 
be coupled with a value of D\ = 9 x 10 8 eV/cm and 
all of the first-order phonons are assumed to be coupled 
with D\\ = 5.6 eV. (This value is consistent with the 
results given in Ref. [57h . The first Born approxima- 
tion result for the total electron-bulk phonon scattering 
rate for p-type silicon with N a = 10 15 cm" 3 , N a = 10 12 
cm" 2 , and T=300 K, with (thick line) and without (thin 
line) the inclusion of the correction term for the first or- 
der process, and for the lowest subband of the unprimcd 
ladder of subbands, is given in Figure [T3"l We see that, 
throughout the whole energy range, there is an increase 
of approximately 10% of the total electron-bulk-phonon 
scattering rate due to the correction term introduced pre- 
viously that could lead to mobility reduction. The same 
trend was also observed for the higher lying subbands. 



3. Polar optical phonon scattering 

Polar optical phonon (POP) scattering is a very 
strong scattering mechanism in polar semiconductors, 
such as GaAs. (It is absent in nonpolar materials, such 
as Si.) Since our focus here is on Si-based structures, we 
will not discuss POP scattering in detail, just note that 
the strength of this mechanism lies in its electrostatic 
nature. Namely, an optical phonon in a polar binary 
semiconductor such as GaAs displaces the two atoms in 
the unit-cell basis with respect to one another, and lead 
to a modification of the dipole moment associated with 
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FIG. 13. Total electron-bulk phonon scattering rate for the 
electrons in the lowest unprimed subband calculated within 
the first Born approximation. 

the unit cell. Overall, the displacement field associated 
with the propagation of an optical phonon gibes rise to 
an electric field, which is what electrons scatter from. 
Detailed derivations of the POP scattering rates can be 
found in many texts££^. Here we will just give the ma- 
trix element squared for POP scattering: 

\(n\U^\m)\^ 2 ^ L + ° ql) {No + \T\) (51) 

x(^-^u)) \ F nm(q z )\ 2 6(k n -k\ { ±q n ) 

where F nm (q z ) is given by Eq. (|45|) , ujlo is the longitudi- 
nal polar optical phonon frequency and and e(0) are 
the high frequency and the low frequency dielectric con- 
stants, respectively. The scattering rate is then found via 
integration over the final states and the final expression 
can be found in Refs. [3 and [33. 

B. Scattering by Electrostatic Interactions 

1. Coulomb Scattering 

Scattering associated with charged Coulomb centers 
near the plane of the 2D electron gas in MOS devices can 
be separated into contributions from the depletion layer, 
the interface charge and the oxide charge. An exten- 
sive discussion of the role of multiple-scattering contribu- 
tions to the electron mobility of a doped semiconductor 
and the apparent difficulties with the impurity averaging 
is presented in the papers by Moored and Kohn and 
Luttinger,^^ The expressions for the potential due to 
a single charge located in the region of interest, in which 
the image term is properly included, are given by Stern 
and Howard. 61 Here, we will discuss it in a way that is 
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suitable for a many-subband treatment. Using the usual 
method of images,— one easily finds that, in the pres- 
ence of dielectric medium, the potential due to a charged 
center located at = (r.;, Zj) equals 



Ui(r,z) = 



1 



4ttk ^(r - r 4 ) 2 + (z - z,) 2 



(52) 



for Zi < 0, and 



1 



^8c 



^ox 



1 



v/(r - r,) 2 + (z - z,) 2 
1 



y/(r - r 4 ) 2 + (z + z 4 ) 2 



(53) 



for Zi > 0. In Eqs. ((52j) and (f53|l. k = 0.5(e sc + e ox ) is the 
average dielectric constant at the interface. The deple- 
tion charge scattering occurs due to the ionized charges 
in the depletion layer. Using Eq. ((53)) . we find that the 
matrix element squared for scattering between subbands 
n and m due to the depletion charge is equal to 



\{n\U d ^(q)\m)\ 2 =\U^\q)\ 2 
^dcpi (^J A 2 nm {q) r o °° dzi Ol m {q, Zi ), 



(54) 



where iVdcpi is the depletion charge density, R; = (rj, Z\) 
is the location of an arbitrary charge center in the de- 
pletion region, and A nm (q) and O nm (q, zi) are the form 
factors due to the finite extension of the electron gas in 
the quantization direction, of the form 



A n m(q) 



dzi/i n (z)e qz i> m {z) 



(55) 



and 



O nm (q, z^ = 0.5 ( J 
+0.5 (l + £ 
respectively, where 



0.5 1 



-qzi 



A nm (q) 



A nm (q) 



,(+) 



dz ip n (z)e ±qz Tjj m (z) 



,(56) 



(57) 



In the above expressions, q is a wavevector in the plane 
parallel to the interface (xy-plane in our case). 

Near the Si/SiC>2 interface, there are always many 
Coulomb centers, due to the disorder and defects in the 
crystalline structure in the neighborhood of the interface. 
They are associated with the dangling bonds and act 
as charge-trapping centers which scatter the free carriers 
through the Coulomb interaction. Using Eq. (|5"2"|) , we 
find that the matrix element squared for the scattering 
from a sheet of charge with charge density located in 
the oxide, at distance Zj (z, < 0) from the interface is 



(n\U ij (q)\m) 



fe(q)l 5 



(58) 









q 



For interface-trap scattering Zj = 0. By similar argu- 
ments, one finds that the matrix element for scattering 
from the oxide charge, with charge density N ox , is 



l(n|[/°*(q)|m)r 



|CC.(q) 



(59) 



= N nx I — 
2k 



2 r 



o,2qd 



2q 



where d ox is the oxide thickness. 

For GaAs/AlGaAs heterostructures, Zj < and re- 
mote Coulomb scattering, in addition to polar optical 
phonon scattering, dominates the low-field electron mo- 
bility. Since there are no charges in the depletion layer 
in this material system (GaAs is intentionally left un- 
doped), direct Coulomb scattering mechanism can safely 
be ignored. In addition to the smaller effective mass 
of electrons in the GaAs system, the absence of direct 
Coulomb scattering is a one of the main reasons for the 
observation of the very high mobility in GaAs modula- 
tion doped heterostructures. 



2. Surface-roughness scattering 

This scattering mechanism is associated with the in- 
tcrfacial disorder and depends upon the oxidation tem- 
perature and ambient as well as post-oxidation anneal 
and removal of the wafer from the furnace. Early theories 
of surface roughness were based on the Boltzmann equa- 
tion in which the surface is incorporated via boundary 
conditions into the electron distribution functionj 63 ' 64 i 65 
The first quantum-mechanical treatment of the problem 
was given by Prange and Nee^S Subsequently, the theory 
followed two different paths. The basic idea of the first 
approach is to incorporate the variations in the confin- 
ing potential of the rough surface as a boundary condi- 
tion on the Hamiltonian of the system. Since there is no 
simple perturbation theory to treat arbitrary changes in 
the boundary conditions, the problem of a free-electron 
Hamiltonian with complicated boundary conditions is 
then transformed by an appropriate coordinate trans- 
formation into a problem with simpler boundary condi- 
tions (i.e. into a problem where we have flat surfaces). 
This coordinate transformation technique has been pro- 
posed by Tesanovic et al^ and was later used by Trivedi 
and Ashcroft.— As a consequence of this transformation, 
the Hamiltonian of the system now has additional terms 
that play the role of potential interaction terms. These 
additional terms are treated by perturbative techniques, 
which are valid when the roughness of the surface is small 
compared to the thickness of the well. 

In the second approach^ the effect of the surface 
roughness is taken into account through a random local 
potential term 

V 6[-z + A(r)] - V 0(-js) ~ U <S(z)A(r), (60) 

which is then treated perturbatively (O is the step func- 
tion). The random function A(r) is a measure of the 
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roughness and is most conveniently expressed in terms 
of the autocovariance function of A(r). The power spec- 
trum S(q) is the two-dimensional Fourier transform of 
the autocovariance function of A(r). For the Gaussian 
correlated roughness that is usually assumed ) 69 ' 70 ! 71 ' 72 ! 73 
the power spectrum is given by 



S G (q) =7rA^exp 



q 2 e 



(61) 



Parameters A and ( characterize the r.m.s. height of 
the bumps on the surface and the roughness correlation 
length, respectively. Goodnick et alJ^ have made an 
extensive analysis of high-resolution transmission elec- 
tron microscopy (HRTEM) measurements to test the as- 
sumption of Gaussian correlation. They found that expo- 
nential correlation describes roughness much better than 
Gaussian correlation irrespective of growth conditions. 
Roughly speaking, it means that the interface may be 
regarded as consisting of terraces of a few nanometers in 
size separated by atomic steps of a few tenths of nanome- 
ters, as shown in Fig. [T?J This result was later confirmed 
by Atomic Force Microscope (AFM) measurements.^ 
The power spectrum for the exponential correlation is 
given by 



Se (q) = 



7rA 2 C 2 



1 



g 2 C 2 

2 



3/2 ' 



(62) 



A generalization of the result given in Eq. (|62p is a self- 
affine roughness correlation function, which in 2D takes 
the form 



SsA(q) 



ttA 2 ( 



2/-2 



n+l : 



(63) 



where n > is an exponent describing the high-g fall-off 
of the distribution. It reduces to exponential correlation 
for n = 0.5. 

For identical roughness parameters, the Gaussian 
spectrum decays slower for small wavevectors and then 
falls to zero rapidly for large wavevectors. The exponen- 
tial model also leads to a rougher interface due to the 
tails in the spectrum, which allows for short-range fluc- 
tuations to be considered as well. For n < 0.5 and small 
values of q, the power spectrum of the self-affine model 
decays faster compared to the previous two models, but 
then falls slowly for large wavevectors. This essentially 
means that, in this regime, it also allows for short-range 
fluctuations to be considered. For large exponents, the 
power spectral density of the self-affine model approaches 
the one for the Gaussian model. In general, the matrix 
element for scattering between subbands n and m for this 
scattering mechanism is of the form 
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FIG. 14. (Top panel) High-resolution transmission electron 
micrograph of the interface between Si and Si02- The oxide 
is in the top half of the picture, while the rows of Si atoms 
can be observed in the bottom half. The image is a lattice 
plane image lying in the (111) plane, while the interface is a 
(100) plane. Reprinted with permission from Ref. [74, S. M. 
Goodnick et al, Phys. Rev. B 32, 8171 (1985). © 1985 The 
American Physical Society. (Bottom panel) Relevant dimen- 
sions for the steps occurring at the interface. 



For large Vq, matrix element T nrn reduces to 
h 2 dip n dip m 



r (0) = 



(65) 



\(n\U sr ( q )\m)\ 2 = S(q)T 2 nm (q). 



(64) 



2m z dz dz 

dz I ip n (z) dV } Z ^ ip m (z) 
o I dz 

- c m ——ii m {z)+t n ip n (z)—— 
dz dz 



The last expression is the result obtained by Prange and 
Nee^ 6 - Matsumoto and Uemura 7 - 6 . calculated that in the 
electronic quantum limit, T nrn = eE av , where E av oc 
4iV s + Ndepi {N s and Nd ep i are the inversion layer and 
the depletion region sheet charge densities, respectively). 

The change in the potential energy of the system due 
to surface-roughness was corrected by Ando^ by con- 
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sidering the change in the electron density distribution 
and the effective dipole moment of the deformed Si-Si02 
surface. The later scattering rate becomes 



r (a) - r (0) + — — 

esc £st 

x <N dep i + N s 



~A nm (q) 



(66) 



where Tnm is given by Eq. (|63|) . Due to the presence of 
the dielectric medium, one needs to correct the expression 
given in Eq. (|66|) with the contribution of the image term 



j^image 
nra 



(?) = 



167re sn e. 



Kx(qz) le s 



-ox Jo 



qz 



2e. 



l K (qz) 



(67) 



i>m(z), 



where K and K x are the modified Bessel functions. An 
additional complication associated with the finite oxide 
thickness, which may further reduce the mobility via 
scattering with remote roughness, will be ignored in the 
present treatment. 



3. Electron- Electron Interaction 

We consider the Coulomb interaction between an elec- 
tron with wave vector k in subband n and a second elec- 
tron with wave vector k2 in subband n-i- The final states 
of these two electrons are k and n for the first electron, 
and k 2 and n 2 for the second electron. The scattering 
rate for this energy-conserving binary collision process 
may be written a o 40 ' 78 



/ 1 1 1 1 

P ee In, k; n 2 , k 2 — > n , k ; n 2 , k 2 



= ^\M^\ 2 6 
n 



£„,(k)+£ n ,(k 2 )-£ n (k)-£ n2 (k 2 ) 

(68) 

where M ee is the matrix element between these two ini- 
tial and final states. With the inclusion of the exchange 
effect for indistinguishable particles, the square of the 
matrix element becomes 7 -^ 



\M e 



+ 



n ,k ;n 2 ,k 2 \H ee \ n, k; n 2 , k 2 

n 2 ,k 2 ;n,k \H ee \ n, k; n 2 , k 2 
n ,k ;n 2 ,k 2 \H ee \ n, k; n 2 , k 2 

n 2 ,k 2 ;n,k |if ee | n, k; n 2 , k 2 



(69) 



n ,k ;n 2 ,k 2 |i? ee | n, k; n 2 , k 2 



|V(q)r 



F e f , (q) 

n n 2 nn-2 v y 



(70) 



where 



V(q) denotes the Fourier transform of the unscreened 
Coulomb potential and q = k — k is the exchanged mo- 
mentum of the first electron with q = |k — k|. The form 
factor is defined by 

f l; M! (q) = J J dzdz 2 Xn'( Z )X n ' 2 (z2)Xn(z)Xn 2 (z 2 ) 

x exp(— q\z — z 2 \). 

This form factor appears also in the random phase 
approximation (RPA) of the dielectric function discussed 
next. 



4- Screening of Coulomb, Surface- Roughness and 
Electron-Electron Scattering 

It is known that Coulomb and surface-roughness scat- 
tering affect significantly the inversion layer electron mo- 
bility, particularly at low and high inversion charge den- 
sities. As we said earlier, the scattering potentials for 
these two dissipative processes are strongly affected by 
the screening of the mobile charges in the inversion layer. 
Therefore, any theory that tries to explain the density 
and temperature dependence of the electron mobility 
must account for these screening corrections. Since the 
calculation of the exact dielectric function of homoge- 
neous electron gas is a formidable problem, various ap- 
proximate solutions for the dielectric function exist in the 
literature^ 1 - 

Some of these have been very successful, because 
of their simplicity (Thomas-Fermi method) or high ac- 
curacy [the random- phase approximation (RPA)]. The 
Thomas-Fermi method is basically the semiclassical limit 
of the Hartree calculation. On the other side, the RPA 
is an exact Hartree calculation of the charge density 
in the presence of the self-consistent field of the ex- 
ternal charge plus electron gas. More precisely, in the 
RPA one includes only the long-range Coulomb interac- 
tion in the dielectric response, leaving out all exchange- 
correlation corrections. It leads to the so-called Lind- 
hard dielectric function that is extensively employed in 
the literature^^^^^^^ 



IV. TRANSPORT IN QUASI-2D SYSTEMS 

The Ensemble Monte Carlo technique has been used 
for over 30 years now as a numerical method to simulate 
nonequilibrium transport in semiconductor materials and 
devices, and has been the subject of numerous books and 
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reviews! 4 ' 89 ! 90 In the application to transport problems, 
a random walk is generated to simulate the stochastic 
motion of particles subject to collision processes in the 
medium. This process of random walk generation may be 
used to evaluate integral equations and is connected to 
the general random sampling technique used in the eval- 
uation of multi-dimensional integrals The basic tech- 
nique is to simulate the free particle motion (referred 
to as free flight) terminated by instantaneous random 
scattering events. The Monte Carlo algorithm consists 
of generating random free-flight times for each particle, 
choosing the type of scattering occurring at the end of 
the free flight, changing the final energy and momentum 
of the particle after scattering, and then repeating the 
procedure for the next free flight. Sampling the parti- 
cle motion at various times throughout the simulation 
allows for statistical estimation of physically interesting 
quantities, such as the single particle distribution func- 
tion, the average drift velocity in the presence of an ap- 
plied electric field, the average energy of the particles, 
etc. By simulating an ensemble of particles representa- 
tive of the physical system of interest, the non-stationary 
time-dependent evolution of the electron and hole distri- 
butions under the influence of a time-dependent driving 
force may be simulated. 

The particle-based picture, in which the particle mo- 
tion is decomposed into free flights terminated by in- 
stantaneous collisions, is basically the same picture un- 
derlying the derivation of the semi-classical Boltzmann 
transport equation (BTE). In fact, it may be shown that 
the one-particle distribution function obtained from the 
random-walk Monte Carlo technique satisfies the BTE 
for a homogeneous system in the long-time limits 



A. Mobility Calculation in Silicon Inversion Layers 
Using 2D Monte Carlo and Including Degeneracy 
Effects 

Multi-particle effects relate to the interaction between 
particles in the system, which is a nonlinear effect when 
viewed in the context of the BTE, due to the depen- 
dence of such effects on the single particle distribution 
function itself. Particle-particle interactions are impor- 
tant in Monte Carlo simulation in establishing or relax- 
ing to an equilibrium distribution function characterized 
by a Maxwell-Boltzman distribution for non-degenerate 
situations, or a Fermi-Dirac distribution when proper ac- 
count for the Pauli exclusion principle is included. Most 
algorithms developed to deal with such effects essentially 
linearize the BTE by using the previous value of the dis- 
tribution function to determine the time evolution of a 
particle over the successive time-step. Multi-carrier ef- 
fects may range from simple consideration of the Pauli ex- 
clusion principle (which depends on the exact occupancy 
of states in the system), to single-particle and collective 
excitations in the system. Inclusion of carrier-carrier in- 
teractions in Monte Carlo simulation has been an active 
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principle 

FIG. 15. Mobility vs. sheet electron density for a silicon in- 
version layer. The wavefunctions of the MOS capacitor struc- 
ture were calculated using SCHRED 39 and fed as input to the 
2D Monte Carlo for the solution of the Boltzmann transport 
equation for confined carriers. 



area of research for quite some time and was discussed 
in conjunction with the description of the scattering 
mechanisms in the previous section. In that discussion, 
we only considered binary collisions but higher multi- 
particle phenomena might be important under high den- 
sity limit, which are usually referred to as the electron- 
plasmon interaction are will not be discussed in this text. 

Another carrier-carrier effect that is of considerable 
importance when estimating, for instance, leakage cur- 
rents in MOSFETs, is impact ionization, which is a pure 
generation process involving three particles (two elec- 
trons and a hole or two holes and an electron) . The Pauli 
exclusion principle requires that the bare scattering rate 
be modified by a factor 1 — /m(k) in the collision integral 
of the BTE, where / m (k) is the one-particle distribution 
function for the state k in band (subband) m after scat- 
tering. Since the net scattering rate including the Pauli 
exclusion principle is always less than the bare scattering 
rate, a self-scattering rejection technique may be used 
in the Monte Carlo simulation as proposed by Bosi and 
Jacoboni 92 for one particle simulation and extended by 
Lugli and Ferry 9 ^ for EMC. In the self-scattering rejec- 
tion algorithm, an additional random number r is gen- 
erated (between and 1), and this number is compared 
to / m (k), the occupancy of the final state (which is also 
between and 1 when properly normalized for the numer- 
ical k-space discretization). If r is greater than / m (k), 
the scattering is accepted and the particle's momentum 
and energy are updated accordingly. If this condition is 
not satisfied, the scattering is rejected, and the process 
is treated as a self-scattering event with no change of 
energy or momentum after scattering. Through this al- 
gorithm, no scattering to this state can occur if the state 
is completely full. 

Using the rejection technique outlined above, we cal- 
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Degenerate statistics 

N s = 2.71x10 13 cm- 2 
E eff = 3.01x1 6 V/cm 



FIG. 16. Distribution functions for a non-degenerate (top) and a degenerate (bottom) quasi-2D electron gas, for given values 
of the sheet electron density and the effective electric field perpendicular to the semiconductor-oxide interface. 



culated the time dependence of the drift velocity along 
a [100] field direction and used a very low value for the 
electric field (lkV/cm) so no heating of the carriers was 
observed. The statistical ensemble consisted of 10,000 
particles that were initially distributed amongst different 
subbands based on the occupancy factors we obtained 
from SCHRED 30 and simulations were performed for 100 
ps. The results of the last 50 ps were then time averaged 
and the mobility was extracted. The variation of the 
carrier drift mobility vs. effective electric field is given 
in Figure [13 We see that our simulator clearly repro- 
duces the experimental data and the Green's function 
results in the mid-to-high effective field region. In these 



simulations, we did not include Coulomb scattering and 
because of that in the low effective field region our simu- 
lation results overshoot the experimental values. It is also 
important to note about a 10% mobility increase in the 
high sheet electron density region (large effective fields) 
because of the inclusion of the Pauli exclusion principle 
which, in turn, leads to closer agreement between the ex- 
perimental data and the Monte Carlo simulations in the 
high effective field region. The distribution functions of 
a non-degenerate and degenerate electron gas are shown 
in Figure [TBI As expected, the degenerate distribution is 
sharper. 



V. QUASI-1D SEMICONDUCTORS 
(NANOWIRES) 



The study of electronic transport in nanowires started 
almost three decades ago with the seminal work by 
Sakaki— in 1980. Sakaki showed that the elastic scat- 
tering in nanowires (NWs) is suppressed drastically due 
to a reduction in the final density of states (DOS) for 
scattering. Electrons in NWs are confined in two trans- 
verse directions and are free to move only along the axis 



of the wire. In the extreme quantum limit (electrons oc- 
cupy only the lowest subband), for an elastic scattering 
process, an electron in the initial state ki can only be 
scattered to the final state kf = —ki accompanied by a 
large change in momentum (2fe,), since the ID constant 
energy surface has only two discrete states. In the case of 
bulk structures, the constant energy surface is a sphere of 
radius ki, so electrons can be scattered to various states 
including those in the vicinity of ki. Sakaki considered 
scattering from remote impurities and showed that the 
mobility as high as 10 s cm 2 /Vs can be achieved in GaAs 
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NWs at low temperatures because of the substantial re- 
duction of the DOS for scattering in NWs. While the 
approximations in this work would underestimate the ef- 
fect of impurities located close to the wire and remote 
impurity scattering is certainly not the only mechanism 
limiting the mobility in NWs, the work of Sakaki pro- 
vided strong impetus for the further study of electronic 
transport in quasi-lD structures. 

In the following year, Arora^ calculated the scatter- 
ing rate due to acoustic phonons and point defects in 
thin rectangular wire under the relaxation time approxi- 
mation (RTA) and found them both to increase with de- 
creasing wire cross section. He also showed that the ratio 
of the conductivity in the wire to the bulk conductivity 
is proportional to the area of the wire cross section. This 
result contradicts the previous result of Sakaki, who as- 
sumed Coulomb scattering from remote impurities alone 
in the calculation of the mobility. Lee and Spector 96 
calculated the impurity-limited mobility by accounting 
for both background and remote impurities using the 
same approximation as in Sakaki's work. They found the 
scattering rate from background impurities to be much 
higher than that from the remote impurities and indepen- 
dent of the wire cross section (probably lost due to the 
delta-function approximation of wavefunctions along the 
transverse directions). They also confirmed the acoustic 
phonon-limited mobility trend in NWs shown by Arora. 

Lee and Vassel^l considered scattering of electrons 
from acoustic phonons, impurities (both remote and 
background), and polar optical phonons and calculated 
the mobility in wires of various cross sections and differ- 
ent temperatures. For the whole temperature range, they 
found the mobility in NWs to decrease with decreasing 
wire cross section. At very low temperatures, where the 
scattering from impurities dominates transport, mobil- 
ity in NWs was found to be higher than that in bulk. 
At room temperature, where phonons dominate electron 
transport, the mobility in wires of cross section smaller 
than 12x12 ran 2 was found to be lower than that of the 
bulk, but for larger cross sections the mobility of the wire 
was greater than in the bulk. 



A. Electronic Bandstructure Modification in 
Nanowires 

All the work discussed above was done on GaAs NWs, 
with the bulk bandstructure, and assuming unconfined 
phonons as in the bulk. Sanders et al.^ in their work on 
the electronic transport in free-standing silicon nanowires 
(SiNWs) , indicated yet another important consequence of 
2D confinement of electrons in NWs - the modification 
of the electronic bandstructure because of the change in 
the dimensionality of the Brillouin zone. The Brillouin 
zone becomes ID in NWs since the crystal structure is 
periodic only along the wire axis. Sanders et al. con- 
sidered SiNWs with axis along [001] and the faces of the 
wire are {110}. The primitive cell from which the wire is 



constructed is shown in Fig. [TTj It contains four silicon 
atoms (as opposed to only two atoms in a bulk silicon 
primitive cell); the length of the primitive cell along the 
axis is a and the transverse dimensions are a/y/2, where 
a is the lattice constant (a = 5.43A). The Brillouin zone 
in SiNWs is ID and it extends from —ir/a to +7r/a, as 
opposed to —2-rr/a to +27r/a in bulk silicon along [100]. 
This is because of the doubling of the length of the unit 
cell along the wire axis. The conduction band in bulk 
silicon is indirect and is composed of six equivalent A 
valleys located at ±0.85 x (2w/a) = ±1.7w/a along each 
of the <100> directions. In case of the SiNWs with axis 
along [001], four of the A4 valleys along the transverse 
directions ([010], [OlO], [100], and [TOO]) are projected 
onto the T point in the ID Brillouin zone and their en- 
ergies are determined by the effective masses along the 
[110] and [TlO] confinement directions. The two A2 val- 
leys along [001] of the bulk Brillouin zone are zone-folded 
to ±0.37r/a in SiNWs and become the off-r states. The 
energy bands derived from these are at higher energies 
than those at the T point since the [001] valleys have a 
lighter effective mass in both confining directions. Thus 
the SiNW becomes a direct bandgap material. 

The work of Sanders et al. paved the way for the 
study of the effects of electronic bandstructure modifica- 
tion on the ballistic transport in silicon nanowire tran- 
sistors (SWNT) by various gro^Zm&MlLmmm 
Using the bandstructure obtained from a tight-binding 
model, Wang et al£2- studied the validity of using the 
bulk effective mass in the calculation of electrical prop- 
erties in SNWTs of different side lengths a. The bulk 
effective mass was found to overestimate both the thresh- 
old voltage (for a < 3 nm) and the ON-current (for a < 5 
ran). Zheng et al*^- calculated the effective masses, val- 
ley splitting between the A4 (at T point) and the A2 (off- 
r point) bands, and the bandgap from the full SiNW 
bandstructure obtained using the tight-binding model. 
All three were found to increase with decreasing wire 
cross section. They also found that the single band effec- 
tive mass equation predicts the bottom of energy bands 
accurately. Bandstructure in SiNWs is also found to de- 
pend on the orientation of the wire axis 1 ^. and the species 
used for surface termination* 1 ^ For all orientations and 
surface terminations studied, SiNWs were found to ex- 
hibit direct bandgap. 

Experimental results on electronic transport in 
SiNWs were first reported by Cui et al. in 2003. The 
mobility in SiNWs was found to be higher than in bulk 
silicon. The increase in mobility was attributed to the 
reduced DOS for scattering as suggested by Sakaki. 94 
In the following year, Koo et al^^ also reported the 
electron mobility in SiNW FETs to be two times higher 
than that in bulk MOSFETs, supporting the results of 
Cui et al. In the same year, Kotlyar et al*^- showed 
that the phonon-limited mobility in SiNWs is much lower 
than that in the bulk silicon because of the increase in 
the electron-phonon wavefunction overlap with decreas- 
ing cross section. Unlike all previous theoretical studies 
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FIG. 17. Crystal structure of an idealized silicon NW. The 
NW unit cell is shown on the right. Each base unit contains 
four atoms. The faces are parallel to the four equivalent {110} 
planes and the wire is oriented along [001]. Reprinted with 
permission from Ref. 193. G. D. Sanders, C. J. Stanton, and 
Y. C. Yang, Phys. Rev. B 48, 11067 (1993). © 1993, The 
American Physical Society. 



of transport in NWs, they considered a detailed multi- 
subband transport and used the electronic wavefunctions 
and subbands obtained from solving coupled Schrodinger 
and Poisson equations self-consistently in the calculation 
of the scattering rates. Wang et al. reported that the 
surface roughness scattering (SRS) becomes less impor- 
tant in the case of ultrasmall SiNWs Jin et al*£ did 
a detailed calculation of both phonon-limited and SRS- 
limited mobility in cylindrical SiNWs and found both to 
decrease with decreasing wire cross section. A similar 
mobility reduction in SiNWs has been reported this year 
by both theoretica l 13 i 108 and experimental groups*^ 



B. Acoustic Phonon Confinement in SiNWs 

Another important consequence of 2D spatial con- 
finement is NWs is the modification of the acoustic 
phonon dispersion in them. In ultrasmall structures 
such as NWs, the acoustic phonon spectrum is modi- 
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FIG. 18. Bandstructures of square SiNWs of three different 
widths. The Brillouin zone is ID with k ranging from —n/a 
to +7r/a. Reprinted with permission from Ref. [9S|, G. D. 
Sanders, C. J. Stanton, and Y. C. Yang, Phys. Rev. B 48, 
11067 (1993). © 1993, The American Physical Society. 



fied due to the mismatch of the sound velocities and 
mass densities between the active layer and the surround- 
ing materialf^i^ in our case Si and Si02- This modifi- 
cation in the acoustic phonon spectrum becomes more 
pronounced as the dimensions of the active layer be- 
come smaller than the phonon mean free path, which is 
around 300 nm in siliconi^ii Pokatilov et a/i 112 i 113 have 
shown that the modification in the acoustic phonon dis- 
persion in nanowires can be characterized by the acoustic 
impedance Q — pV s , where p and V s are the mass density 
and sound velocity in the material, respectively. By con- 
sidering materials with different £, Pokatilov et al. have 
shown that the acoustic phonon group velocity in the ac- 
tive layer is reduced when an acoustically soft (smaller £) 
material surrounds an active layer made of acoustically 
hard (higher £) material. Since Si is acoustically harder 
than Si02, the acoustic phonon group velocity in SiNWs 
with Si02 barriers decreases and results in an increased 
acoustic phonon scattering rate [see Eq. (|7"3[)] . 

The first step in accounting for the acoustic phonon 
confinement is to calculate the modified acoustic phonon 
dispersion. Using the adiabatic bond charge model 115 
(microscopic calculation, accurate but computationally 
involved), Hepplestone and Srivastava 11 - have shown the 
validity of the elastic continuum model (macroscopic cal- 
culation, less accurate but easier to implement) for wire 
dimensions greater than 2.5 nm. Hence, in this work we 
have used the elastic continuum model to calculate the 
modified phonon spectrum. Most of the previous stud- 
ies of acoustic phonon confinement in nanowires have 
used approximate hybrid modes proposed by Morseiii 
(valid for wires with thickness much smaller than the 
width) to calculate the dispersion spectrum. Nishiguchi 
et alr^^ calculated the dispersion spectrum using the xyz 
algorithm 118 and found that the Morse formalism is valid 
only for the lowest phonon subband. Since one acoustic 
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FIG. 19. Confined acoustic phonon dispersion (dilatational 
mode) calculated using the xyz algorithm^ for an 8 x 8 
nm 2 SiNW. Only the lowest 10 phononic subbands are shown. 
Dispersion in the first one third of the first Brillouin zone is 
shown for clarity. Reprinted with permission from Ref. Il3l 
E. B. Ramayya et al., J. Appl. Phys. 104, 063711 (2008). © 
2008, American Institute of Physics. 



phonon subband is surely not enough to accurately de- 
scribe scattering with electrons, in this work, we have 
used Nishiguchi et al.'s approach to calculate the acous- 
tic phonon dispersion, although it is computationally in- 
tensive. The basis functions used to expand the phonon 
mode displacements in Nishiguchi et a/.'s approach are 
powers of Cartesian coordinates, and the number of basis 
functions required to fully describe the modes depends on 
the number of modes required. For an 8 x 8 nm 2 SiNW, 
we have found that the lowest 35 phononic subbands are 
enough to calculate the scattering rate. Also, we found 
that about 176 basis functions are sufficient to fully de- 
scribe the displacement of those 35 phononic modes. The 
number of phononic bands required decreases with the 
decrease in wire cross section. 

Two types of boundary conditions are often used to 
calculate the acoustic phonon spectrum in nanostruc- 
tures: a) The free-standing boundary condition (FSBC) 
assumes that all the surfaces are free, so normal com- 
ponents of the stress tensor vanish at the surfaces, and 
b) The clamped-surface boundary condition (CSBC) as- 
sumes that the surfaces are rigidly fixed, so the displace- 
ment of phonon modes is zero at the surfaces. Generally, 
the CSBC (FSBC) results in higher (lower) phonon group 
velocity than the bulk case^iS For the wires considered 
in this work, neither of these boundary conditions holds 
exactly, since these wires are actually embedded in the 
Si02- Ideally, one needs to solve the elastic continuum 
equation, taking into account the continuity of displace- 
ment and stress at all Si-Si02 interfaces, and then apply 
the boundary conditions at the outer surfaces. But, this 
is almost numerically impossible for the structure con- 



sidered because it would be equivalent to solving the 2D 
Schrodinger equation in a device with the cross section of 
about 800 x 400 nm 2 , three (five if the metal-Si interface 
is included) interfaces along the depth, and two interfaces 
along the width. Donetti et al., in their work on a SiOv 
Si-Si02 sandwich structure, considered continuity of the 
displacement and stress at the interfaces to calculate the 
phonon dispersion*^ They found it to be close to the 
results from FSBC. Therfore, in this work, we have used 
the FSBC to calculate the acoustic phonon spectrum of 
SiNWs. Fig. Q33 shows the calculated acoustic phonon 
dispersion of the lowest 10 dilatational modes for an 8 
x 8 nm 2 SiNW. Apart from these dilatational modes, 
depending on the rotational symmetry of the confined 
acoustic phonon displacement, there are two sets of flex- 
ural modes and one set of torsional modes in SiNWs. A 
detailed description of the s ymm etry of all these phonon 
modes can be found in Ref. Ill4 



VI. SCATTERING IN QUASI-1D ELECTRON 
SYSTEMS 

A. Scattering due to Bulk Acoustic Phonons, 
Intervalley Phonons, and Surface Roughness 

Phonon scattering and the SRS are the most impor- 
tant scattering mechanisms in SiNWs. The SRS was 
modeled using Ando's model^ intervalley scattering was 
calculated using bulk phonon approximation, and the in- 
travalley acoustic phonons were treated in both the bulk- 
mode and confined-mode approximations. Since the wire 
is very lightly doped, the effect of impurity scattering was 
not included. Nonparabolic band model for silicon, with 
the nonparabolicity factor a — O.beV^ 1 , was used in the 
calculation of scattering rates. A detailed derivation of 
the ID scattering rates is given in Appendices of Ref. 
IT3I . Here, for brevity, only the final expressions for the 
scattering rates are given. 

For an electron with an initial lateral wavevector k x 
and parabolic kinetic energy £fc x = ?t 2 fc 2 /(2m*) in sub- 
band n [with subband energy £ n and electron wavefunc- 
tion ip n (y, z)], scattered to subband m [with subband en- 
ergy £ m and electron wavefunction il) m {y,z)], the final 
kinetic energy Ef is given by 

x /l + 4a£ fc - 1 

8j=£ n -8 m + ^ — ^ + ?uj, (72) 

za 

where Tilj = for elastic (bulk intravalley acoustic 
phonon and surface roughness) scattering, Tiuj = iTiwo 
for the absorption/emission of an approximately disper- 
sionless intervalley phonon of energy fajQ, while in the 
case of confined acoustic phonons (below) the full phonon 
subband dispersion is incorporated. 

The intravalley acoustic phonon scattering rate due 



21 



to bulk acoustic phonons is given by 
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(k x ) 
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(1 + 2a£ f ) 

y/£f(l+<*£ f ) 



where s a c is the acoustic deformation potential, p is the 
crystal density, V s is the sound velocity, and O is the 
Heaviside step-function. T> nm represents the overlap inte- 
gral associated with the electron-phonon interaction (the 
so-called electron-phonon wavefunction integral, W 6 ) and 
is given by 



\ip n {y,z)\ 2 \ip m (y,z)\ 2 dydz. (74) 



For scattering mediated by short wavelength acoustic 
and optical phonons, the intervalley phonon scattering 
rate is given by 
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(75) 



where S„ is the intervalley deformation potential, and 
V nm is defined in (|74|) . The approximation of dispersion- 
less bulk phonons of energy Tlloq was adopted to describe 
an average phonon with wavevector near the edge of the 
Brillouin zone and Nq = [exp(huJo / k B T) — 1] 1 is their 
average number at temperature T. 

Assuming exponentially correlated surface 
roughnessi 2 ^ and incorporating the electron wave- 
function deformation due to the interface roughness 
using Ando's model,— the unscreened SRS rate is given 
by 



A 2 A 
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(76) 



Q(£ f ), 



where A and A are the r.m.s. height and the correlation 
length of the fluctuations at the Si-Si02 interface, respec- 
tively. q x — k x ±k' x is the difference between the initial 
(k x ) and the final (k' x ) electron wavevectors and the top 
(bottom) sign is for backward (forward) scattering. The 
SRS overlap integral in Eq. (f75|) due to the top interface 
for a silicon body thickness of t v is given by 
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The SRS overlap integral was derived assuming the in- 
terfaces to be uncorrelated. For the bottom interface, 



the integration should be performed from the bottom in- 
terface to the top interface and the integral for the side 
interfaces can be obtained by interchanging y and z in 
Eq. (J77]). The first term in Eq. (JTTJ) is the confinement- 
induced part of the SRS and it increases with decreasing 
wire cross section. This term does not depend on the 
position of the electrons in the channel and hence re- 
sults in high SRS even at low transverse fields from the 
gate. The second and third terms in Eq. (1771) depend on 
the average distance of electrons from the interface, so 
they contribute to the SRS only at high transverse fields 
from the gate. For wires of cross section smaller than 
5x5 nm 2 , major contribution to SRS comes from the 
confinement-induced term in Eq. ([77| . and it increases 
rapidly with decreasing wire cross section. 

Scattering rates given by Eqs. Q73[) - f76"l) are cal- 
culated using the wavefunctions and potential obtained 
from the self-consistent Poisson-Schrodinger solver. Pa- 
rameters used for calculating the intervalley scattering 
were taken from Ref. 1 12 ll the acoustic deformation po- 
tential was taken from Ref. [13, and A = 0.3 nm and A = 
2.5 nm were used to characterize the SRS due to each of 
the four interfaces. The SRS parameters were obtained 
by fitting the mobility of an 8 x 30 nm 2 SiNW in the 
high transverse field region (where the SRS dominates) 
with the corresponding mobility observed in ultra-thin 
SOI of similar thicknessJ^ 2 - 



B. Scattering due to Confined Acoustic Phonons 

The modification of the acoustic phonon dispersion 
due to confinement, shown in Fig. 1191 implies that 
the linear dispersion and elastic scattering approxima- 
tion can no longer be used in calculating the scattering 
rate. The modified scattering rate which takes into ac- 
count confined acoustic phonon modes is given by 
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where q x is the lateral wavevector of the acoustic phonon, 
g(q x ) = £ — £ / ^fhuij(q x ), q Xl and q X2 are the two possible 
roots of g{q x ) = 0, and g'(q Xl ) and g'(q X2 ) are the deriva- 
tives of g{q x ) with respect to q x evaluated at q Xl and 
q X2 , respectively. Index J stands for the different acous- 
tic phonon modes and Nj qa . is the number of acoustic 
phonons of energy Tiujj(q x ). Overlap integral C nm (J,q x ) 
and the total energy of the electron before (£) and after 
{£') scattering are defined in the appendix of Ref. [l3|. 
Xjq w . is the eigenvector corresponding to ojj(q Xi ) (for 
details, see Appendix of Ref. [HI). 

Fig. [201 shows the electron-acoustic phonon scatter- 
ing rate calculated using both bulk-mode and confined- 
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mode approximations. When calculating the electron- 
bulk acoustic phonon scattering rates, acoustic phonon 
dispersion is assumed to be linear. This assumption un- 
derestimates the electron-acoustic phonon scattering rate 
with respect to the calculation with confined phonons be- 
cause the bulk phonon velocity is higher than the veloc- 
ity of confined phonons and the electron-acoustic phonon 
scattering rate is inversely proportional to the acoustic 
phonon group velocity. On average, the confined acous- 
tic phonon scattering rate is about two times the acous- 
tic scattering rate calculated using bulk phonons. Also, 
each of the bulk-phonon intersubband spikes is split into 
two groups of spikes corresponding to absorption and 
emission of different confined phonon modes. In order 
to account for both intrasuband and intersubband tran- 
sitions, all four acoustic phonon modes are included in 
the calculation of the confined acoustic phonon-electron 
scattering rates. But, it should be noted that the domi- 
nant contribution to the scattering rate comes from the 
dilatational modes (for all other modes the overlap in- 
tegral C nm (J,q x ) in Eq. (|78|) vanishes for intrasubband 
scattering due to symmetry). 
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FIG. 20. Electron-acoustic phonon scattering rate of an 8 x 
8 nm 2 SiNW at the channel sheet density of jV 3 = 8.1 x 10 11 
cm -2 , calculated assuming the bulk and confined phonon ap- 
proximations. The electron-bulk acoustic phonon intersub- 
band spikes are at around 20 meV, 52 meV, 85 meV, and 
95 meV and they correspond to the electron scattering from 
the lowest subband to the 2 nd , 3 rd , 4 th , and the 5 th low- 
est subbands. The set of spikes indicated by Ab and Em 
correspond to the absorption and emission of phonons from 
different phonon subbands. Reprinted with permission from 
Ref. GJ, E. B. Ramayya et al, J. Appl. Phys. 104, 063711 
(2008). © 2008, American Institute of Physics. 



VII. TRANSPORT IN QUASI-1D ELECTRON 
SYSTEMS 

In the following, we present recent 
calculation o 13 i 123 i 124 of the electron mobility in rectan- 
gular (Sec. IVTTAl) and square (Sec. \VTlB\i SiNWs, in 
which transport is governed by scattering of electrons 
with acoustic phonons (confined and bulk), intervalley 
phonons, and imperfections at the Si-Si02 interface. We 
consider wires much longer than the electron mean-free 
path under a low lateral electric field. The nanowires 
are formed on thin silicon-on-insulator (SOI), with the 
gate oxide, burried oxide, and bottom silicon substrate 
thicknesses equal to 25 nm, 80 nm, and 700 nm, 
respectively. All gated wires (rectangular and square) 
have 200 nm of oxide on the side. The channel is lightly 
doped to 3 x 10 15 cm -3 and the wires are assumed 
homogeneous and infinitely long. In Sec. IVH C| we also 
briefly discuss transport in highly doped ungated wires 
used for thermoelectric applications, which have only 2 
nm of native oxide all around. 

Bulk-silicon effective mass parameters are used in 
the calculation of the scattering rates. The confined 
acoustic phonon spectrum is obtained by using the xyz 
algorith m 114 ' 118 and the unscreened SRS is modeled us- 
ing modified Ando's model^ that accounts for the finite 
thickness of the silicon layer. The simulator used to cal- 
culate the electron mobility has two components: the first 
is a self-consistent 2D Poisson-2D Schrddinger solver and 
the second is a Monte-Carlo transport kerne L 13 ' 123 ' 124 
The former is used to calculate the electronic states and 
the self-consistent potential distribution across the wire 
and the latter simulates transport along the wire axis. 
The finite barrier at the Si-Si02 interface results in the 
electron wavefunction penetration through the interface 
and into the oxide. The wavefunction penetration is ac- 
counted for by including a few mesh points in the ox- 
ide while solving the Schrddinger equation. ARPACK 
package 1 ^ was used to solve the 2D Schrddinger equa- 
tion and successive over-relaxation (SOR) method was 
used to solve the 2D Poisson equation. The convergence 
of the coupled Schrodinger-Poisson solver is found to be 
faster when the Poisson equation is solved using the SOR 
method than when it is solved using incomplete LU (ILU) 
method. 

A Monte Carlo transport kernel is used to simulate 
electron transport along the axis of the wire under the 
influence of the confining potential in the transverse di- 
rections and a very small lateral electric field along the 
channel. The long wire approximation implies that the 
transport is diffusive (the length exceeds the carrier mean 
free path), and therefore justifies the use of the Monte 
Carlo method 4 -' 1 ^ to simulate electron transport. Elec- 
trons are initialized such that their average kinetic en- 
ergy is (l/2)ksT (thermal energy for ID) and are dis- 
tributed among different subbands obtained from the 
Poisson-Schrodinger solver in accordance with the equi- 
librium distribution. Since the electrons are confined in 
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FIG. 21. Flowchart of the simulator developed to calculate 
the electron mobility in SiNWs. 
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FIG. 22. Variation of the field-dependent mobility with vary- 
ing SiNW width. The wire thickness i s kep t constant at 8 nm. 
Reprinted with permission from Ref. 1 1231 . E. B. Ramayya el 
al, IEEE Trans. Nanotechnol. 6, 113 (2007). © 2007, IEEE. 



two transverse directions, they are only scattered either 
forward or backward; consequently, just the carrier mo- 
mentum along the length of the wire needs to be updated 
after each scattering event. Mobility is calculated from 
the ensemble average of the electron velocities. 4 



A. Gated Rectangular Nanowires 

The variation of the field-dependent mobility with de- 
creasing channel width was investigated on a series of 
SiNWs, while keeping the channel thickness at 8 nm and 
using the bulk phonon approximation. Fig. [55] shows 
the mobility for SiNWs with the widths of 30 nm, 16 nm 
and 8 nm. Two important results regarding the mobility 
behavior in the width range considered can be deduced 
from Fig. [55J (i) the mobility at high transverse fields, 



which is dominated by the SRS, increases with decreas- 
ing wire width and (ii) the mobility at low-to-moderate 
transverse fields, determined by phonon scattering, de- 
creases with decreasing wire width. In the calculation of 
the SRS overlap integral, initially infinite wire thickness 
was assumed as in the original Ando's derivation. With 
t s — > co assumption the integral becomes 



Tnm = I I dydz[ip n (y,z)e y {y,z)^ m (y,z) (79) 
'£ m -£ n \ dtp m (y,z) 



dy 



Phonon scattering variation with decreasing wire width 
is determined by the interplay of two opposing factors: 
(i) reduction of the final density of states for the elec- 
trons to scatter to, and (ii) an increase in the electron- 
phonon wavefunction overlap (|74|) ■ The former results in 
an enhanced mobility, while the latter results in mobility 
degradation. The overlap integral ([74")) increases with a 
decrease in the wire width due to an increase in the elec- 
tron confinement In narrow wires, the increase in the 
electron-phonon wavefunction overlap dominates over the 
density-of-states reduction, resulting in a net decrease 
in the electron mobility at low-to-moderate transverse 
fields. 

The SRS overlap integral given by (177)) has three 
terms. Two are due to deformation of the wavefunction, 
and the third, dominant term depends on the strength 
of the field perpendicular to the interface. Since the field 
normal to the side interfaces is very weak, SRS due to 
these interfaces is much less efficient than scattering due 
to the top and bottom ones. The decrease in SRS with 
decreasing wire width can be understood by following the 
behavior of the average distance (|8D]) of the carriers from 
the top interface: 



(y) = — 

W Ni 



\^i(y,z)\ 2 y dydz, 



(80) 



where Ni is the total line density and N,' v is the line 



density in the i subband of the v valley. In Fig. 
1531 we can see that the carriers are moving away from 
the top interface as the width of the wire is decreased, 
and are therefore not strongly influenced by the inter- 
face. This behavior is also observed in the case of ul- 
trathin double-gate SOI FETs, and is due to the onset 
of volume inversion ^ 27 ! 128 As the SiNW approaches the 
volume inversion limit, carriers cease to be confined to 
the interfaces, but are distributed throughout the silicon 
volume. Fig. [Ml shows the distribution of carriers in a 30 
nm wire and an 8 nm wire at the same effective field; we 
can clearly see that the carriers are confined extremely 
close to the interface in the 30 wide nm wire, whereas 
they are distributed throughout the silicon layer in the 
case of the 8 nm wire. 
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FIG. 23. Variation of the average distance of carriers from 
the top interface (below the gate) for various wire widths as 
a function of the eff ectiv e transverse field. Reprinted with 
permission from Ref. 1 1231 . E. B. Ramayya et al, IEEE Trans. 
Nanotechnol. 6, 113 (2007). © 2007, IEEE. 
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FIG. 24. Electron distribution across the nanowire, for the 
wire width of 30 nm (left panel) and 8 nm (right panel). In 
both panels, the transverse field is 1 MV/cm, the wire thick- 
ness is 8 nm, and the color scal e is in xl0 19 cm -3 . Reprinted 
with permission from Ref. 1 1231 . E. B. Ramayya et al., IEEE 
Trans. Nanotechnol. 6, 113 (2007). © 2007, IEEE. 



B. Gated Square Nanowires 

In this section, we emphasize the importance of acous- 
tic phonon confinement in SiNWs, and vary the cross 
section of the wire to investigate the effect of increasing 
spatial confinement on electron mobility. The SRS over- 
lap integral given by Eq. (|77j) is used to calculate the 
SRS to account for the finite thickness and width of the 
wire. 

The electron mobility in an 8 x 8 nm 2 wire, with and 



without phonon confinement, is shown in Fig. 1251 In 
the low transverse field regime, the mobility calculated 
with confined acoustic phonons is about 10 % smaller 
than that obtained with bulk phonons. This clearly indi- 
cates that confined acoustic phonons need to be properly 
included in the study of electrical transport in SiNWs. 
The mobility values for 8 x 8 nm 2 are very close to the 
experimentally observed mobility in an ultra-thin SOI of 
similar thickness*^ In the rest of the section, acoustic 
phonons are treated under confined mode approximation. 
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FIG. 25. Variation of the field- dependent mobility for an 8 
x 8 nm 2 SiNW assuming bulk acoustic phonons (solid line) 
and confined acoustic phonons (dashed line). Reprinted with 
permission from Ref. 1 13 . E. B. Ramayya et al., J. Appl. Phys. 
104, 063711 (2008). © 2008, American Institute of Physics. 
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FIG. 26. Variation of the electron mobility with SiNW cross 
section at three different transverse fields. Reprinted with 
permission from Ref. [l3|, E. B. Ramayya et al., J. Appl. Phys. 
104, 063711 (2008). © 2008, American Institute of Physics. 

To determine the cross sectional dependence of elec- 
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tron mobility in SiNWs and to understand the confine- 
ment effects on the spatial and fc-space distribution of 
electrons, the cross section of the wire was varied from 
8x8 nm 2 to 3 x 3 ran 2 . The variation of the electron 
mobility with decreasing wire cross section at a low (1.4 
xl(T 2 MV/cm), moderate (2.4 xlCT 1 MV/cm) and high 
(1.04 MV/cm) transverse field is plotted in Fig. [26] 
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FIG. 27. Variation of the electron population at a low trans- 
verse field (1.4 xlO -2 MV/cm). The solid line shows the total 
population of the electrons in the lowest subband in each of 
the four A4 valley pairs and the dashed line shows popula- 
tion of the A2 valley pair, with varying spatial confinement. 
Reprinted with permission from Ref. fl3l . E. B. Ramayya et 
al, J. Appl. Phys. 104, 063711 (2008). © 2008, American 
Institute of Physics. 



E 
c 



Q_ 
CD 
"D 

CD 31 

£ 

O) 

c 
o 

CO 2 
o 4 

"55 

b 



OlOlO 



□ 



1 



Distance along the width [nm] 



FIG. 28. Electron density across the nanowire at a high 
transverse field (1.04 MV/cm). When the cross section is 
reduced from 8x8 nm 2 (bottom right panel) to 3 x 3 nm 2 
(top left), the onset of volume inversion is evident. The color 
scale is in 5 xlO 18 cm" 3 . Reprinted with permission from 
Ref. GJ, E. B. Ramayya et al., J. Appl. Phys. 104, 063711 
(2008). © 2008, American Institute of Physics. 



1. Subband Modulation 

One of the important factors that determine the en- 
ergy and occupation probability of a subband in each 
of the A@ valleys (equivalent in bulk silicon) is the ef- 
fective mass in the direction of confinement. For the 
SiNWs considered, the confinement is along the y and z 
directions and the electrons are allowed to move freely 
in the x direction; consequently, the conductivity effec- 
tive mass for the valley pairs with minima on the x, y, 
and z axes are mi, rat, and m t , respectively, while their 
subband energies are proportional to l/yVi 2 , l/y/m t mi, 
and 1/ y/mtmi, respectively. Since m t < mi, the sub- 
bands in the valley pair along x are higher in energy 
than those in the valley pairs along y and z. So the sub- 
bands split into those originating from the twofold de- 
generate A2 (the valley pair along x) and those originat- 
ing from the fourfold degenerate A4 valleys (the valley 
pairs along y and z). Upon increasing spatial confine- 
ment by decreasing the wire cross section, the subbands 
in different valleys are pushed higher up in energy, and 
consequently only a few of the lowest subbands in each 
of the valley pairs get populated with electrons. Fig. [27] 
shows the depopulation of the higher A 2 valley subbands 
with increasing spatial confinement: since the lowest sub- 
bands in the A4 valleys are lower in energy than those in 
the A2 valleys, under extreme confinement A2 subbands 
get completely depopulated, and only the lowest A4 sub- 
bands are populated. Splitting of the valley degeneracy 
and modification of the subband energies in different val- 
ley pairs due to spatial confinement, followed by depop- 
ulation of the higher subbands, are together termed sub- 
band modulation*^ Although subband modulation re- 
sults in enhanced electron mobility in ultrathin-body SOI 
MOSFETs, 122 the rapid increase in SRS for wire cross 
section below 5x5 nm 2 suppresses the effect of sub- 
band modulation in ultra-small SiNWs. In our previous 
workfi2i we did observe a small enhancement in mobility 
for wires of cross section around 4x4 nm 2 , but in that 
study we did not consider the confinement induced term 
in SRS. This has been shown to be the dominant term in 
determining the SRS in ultra small SiNWs 



2. Volume Inversion 

As the cross section of the SiNW decreases, the chan- 
nel electrons are distributed throughout the silicon vol- 
ume as opposed to just within a thin channel at the Si- 
Si02 interface in conventional MOSFETs. The transition 
from surface inversion to volume inversion occurs grad- 
ually and the cross section at which the entire silicon is 
inverted depends on the electron sheet density*^ Fig. 
|2"81 shows the variation of the electron density across the 
wire with varying wire dimensions. When the cross sec- 
tion is decreased from 8x8 nm 2 to about 6x6 nm 2 , 
the onset of volume inversion results in an increase in 
the average distance of the electrons from the interfaces. 
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But when the wire cross section is around 6x6 nm 2 , 
the silicon layer is fully inverted, so further decrease in 
the cross section simply results in a decrease of the aver- 
age distance of the electrons from the interfaces, thereby 
resulting in more SRS. Consequently, for wires with the 
cross section smaller than 6x6 nm 2 , volume inversion 
does not offer an advantage to electronic transport. 



3. Mobility Variation with the SiNW Cross Section 

Fig. [26] shows the mobility variation when the SiNW 
cross section is varied from 8x8 nm 2 to 3 x 3 nm 2 . 
The three curves correspond to three different transverse 
effective fields. Irrespective of the effective field from the 
gate, the mobility decreases with increasing spatial con- 
finement mainly due to the monotonic increase in the 
confinement-induced part of SRS and the intrasubband 
phonon scattering. At low and moderate effective fields 
from the gate, the mobility is determined by the scatter- 
ing from phonons and the confinement-induced part of 
SRS, whereas at high fields, mobility is limited by SRS 
(all three terms in the SRS overlap integral, given by 
Eq. ([77)) play crucial roles). 

At high fields from the gate, when the wire cross sec- 
tion is reduced from 8x8 nm 2 to 5 x 5 nm 2 , the first 
term in the SRS overlap integral increases, whereas the 
second and the third terms decrease due to the onset of 
volume inversion (Fig. I28p . Consequently, the mobility 
shows a very small change for these cross sections. But, 
when the wire cross section is smaller than 5x5 nm 2 , 
all the terms in the SRS overlap integral increase with 
decreasing wire cross section resulting in a monotonic 
decrease in mobility with increasing spatial confinement. 

At low and moderate fields from the gate, with de- 
creasing wire cross section: intrasubband phonon scatter- 
ing increases due to the increase in the electron-phonon 
overlap integral; intersubband scattering and intervalley 
phonon scattering decreases due to subband modulation; 
SRS increases due to the increase in the first term in the 
SRS overlap integral with increasing confinement. Over- 
all the mobility decreases with decreasing wire cross sec- 
tion. It should be noted that at low transverse fields, the 
mobility variation for wires of cross section larger than 7 
x 7 nm 2 becomes very small. This is because of the in- 
terplay between a simultaneous increase in intersubband 
scattering (number of occupied subbands increases) and 
a decrease in intrasubband scattering (electron-phonon 
overlap integral decreases) with increasing wire cross sec- 
tion. A similar weak dependence of the electron mobility 
on the wire cross-sectional dimension has been reported 
by Jin et al~2. for cylindrical SiNW of diameter above 6 
nm. 
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FIG. 29. Variation of the electron mobility with the wire 
cross section in ungated square SiNWs. The trend in the 
mobility variation is similar to the trend observed in gated 
wires at a very low transverse field. 



C. Ungated Silicon Nanowires 

In contrast to the wires considered so far, SiNWs used 
for thermoelectric applications are ungated. They are 
also very highly doped to increase the electrical conduc- 
tivity. In obtaining the self-consistent wavefunctions and 
potential in ungated wires, the electric field and the wave- 
functions are forced to zero at the air-Si02 interfaces. 
The silicon channel is assumed to be doped to 1.6 x 10 19 
cm~ 3 n-type with arsenic, and the native oxide thickness 
is 2 nm all around. Due to the high doping concentration, 
scattering from ionized impurities is expected to play a 
crucial role in determining the mobility, so in addition to 
the SRS and phonon scattering considered before, scat- 
tering due to impurities is also included in the mobility 
calculation. 

The electron-impurity scattering rate is given by 

v im P(v , gVjWm (l + 2qg / ) 
nm 1 x) 16YW4 y/e f (l + a£ f ) 

x J dRI 2 m (q±,R), (81) 

where Z — 1 is the number of electrons donated by the 
impurity atom, N a is the doping density, e S i is the di- 
electric constant of silicon, e is the electron charge, m 
is the electron mass, a is the non-parabolicity factor, 
£f is the final kinetic energy as defined in Eq. (|T2"j) . 
Qx = k x i k x is the difference between the initial and 
final electron wavevectors, and R is the position of the 
impurity in the wire cross section. The impurity overlap 
integral I 2 m (<g , R) is defined as 
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l nm (q,R) = J ip n {y, z)K (q,r,R)tp m (y, z)dydz, (82) 
where 

V(r- R) 2 + x a 

Ko(«,r,R)= / dx. (83) 

7 >/(r - R) 2 + x 2 

Since the doping level is in the degenerate limit, the Pauli 
exclusion principle is included in the Monte Carlo simu- 
lation via the rejection technique described earlier (Sec. 
IIV[) : an electron is allowed to undergo a given scatter- 
ing event only if the final state is empty, and if not, the 
scattering event is treated as self-scattering.^^ Electron 
mobility in ungated wires is found to decrease with de- 
creasing wire cross section (Fig. 129]) . as previously seen 
in gated wires at low transverse fields (Fig. 121))) . 

VIII. CONCLUDING REMARKS 

In this paper, we have reviewed the semiclassical 
transport description in Q2DEGs and QlDEGs. Spatial 
confinement in Q2DEGs and QlDEGs leads to forma- 
tion of low-dimensional subbands, and a system of cou- 
pled Schrodinger and Poisson equations in the effective 
mass approximation needs to be solved to obtain the sub- 
band energies. Subbands are then populated according 
to the Boltzmann transport equation, whose state-of-the- 
art numerical solution is based on the ensemble Monte 
Carlo technique. Ensemble Monte Carlo remains fairly 
robust for room-temperature transport description down 
to deep-submicron scales. 

Representative examples - the silicon MOSFET and 
silicon nanowires with varying cross section - were ana- 
lyzed in detail. As the MOSFET oxide thickness is scaled 
to below 10 nm, quantum confinement of inversion charge 
leads to an appreciable inversion layer capacitance in se- 
ries with the oxide, so the total gate capacitance is low- 
ered. The influence of the inversion layer capacitance on 



the threshold voltage shift and device transconductance 
has been reliably modeled by self-consistent Schrodinger- 
Poisson solvers. Furthermore, the field-dependent mobil- 
ity curve for the silicon MOSFET obtained in experiment 
has been reproduced very well by using a Schrodinger- 
Poisson-Monte Carlo transport simulation of the inver- 
sion layer Q2DEG. The importance of accounting for the 
degeneracy in the carrier statistics was clearly demon- 
strated. 

In ultrathin SiNWs, both electrons and acoustic 
phonons experience 2D confinement. In wires surrounded 
by Si02, an acoustically softer material, the acoustic 
phonon group velocity is lowered to almost half of its bulk 
silicon value, and leads to enhanced electron-acoustic 
phonon scattering rates. The electron mobility calcu- 
lated while accounting for the modification to the acous- 
tic phonon spectrum due to confinement is about 10% 
lower than the mobility calculated with bulk acoustic 
phonons. For very thin wires (below the 5x5 nm 2 cross 
section), the mobility decreases monotonically with in- 
creasing spatial confinement and becomes virtually in- 
dependent of the transverse electric field. This occurs 
primarily due to the increase in the field-independent, 
confinement-induced part of the SRS, and second due 
to the increase in intrasubband phonon scattering. In 
contrast to bulk MOSFETs, in which the SRS plays an 
important role only for high fields from the gate, elec- 
trons in very thin SiNWs are strongly influenced by the 
roughness regardless of the transverse field. This find- 
ing is important for field-effect transistors with multiple 
gates, as well as for ungated ultrathin wires used for ther- 
moelectric applications or interconnects. 
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